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Análise teorica das caracteristicas mecânicas 


de um escoamento radial numa turbomáquina centrifuga 


RESUMO 


Neste trabalho são apresentadas, de uma forma 
condensada, as deduções das equações gerais das velo- 
cidades e acelerações e energia num escoamento sem 
perdas no interior da roda de uma turbomáguina radial. 
É, seguidamente. feita a aplicação das equações estabe- 
lecidas a vários escoamentos: escoamento a energia 
constante e escoamentos com aumento de energia (velo- 
cidade relativa constante, inclinação redial constante e pá 
plana) e, finalmente, ssccamentos com entrada e saida 
impostas: pá circular e variação linear da inclinação 
tangencial. 


1 — INTRODUÇÃO 


E com o objectivo de analizar diferentes hipóteses 
da escoamentos em turbomáquinas radiais que o presente 
trabalho é desenvolvido, procurando avaliar quantitativa- 
mente os valores das velocidades e acelerações em jogo, 
variação de pressão ao longo da pá e energia total do 
escoamento. 

Os escoamentos estudados são: escoamentos a 
energia constante como caso limite, escoamento se- 
gundo uma espiral logarítmica, escoamento com veloci- 
dado relativa constante e, pá plana. São também es- 
tudados os escoamentos em torno de uma pá circular 
8 escoamento com variação linear da velocidade relativa. 
O estudo é ilustrado com uma aplicação numérica. 

A consideração de escoamentos unidimensionais é, 
como se sabe, demasiado simplista mas nem por isso é 
menos importante o seu estudo pois que é a partir deles 
a aplicando certas correcções que se prevê o desenvol- 


ANTONIO A. 5. DE BARROS 


SUMMARY 


This study establish the velouty and aceleration 
equations for the flow in a centrifugal turbomachine, 
wich are applied to some types of flow: without and 
with energy increase. 


vimento dos escoamentos nas turbomáquinas radiais. 

As relações apresentadas visam, essencialmente, 
uma aplicação directa a perfis de pás bem conhecidos 
que poderão ser assim convenientemente comparados. 

Todos os conceitos, embora muito simples, são 
sempre abordados com a preocupação de evidenciar a 
aplicação dos teoremas da Mecânica aos escoamentos 
unidimensionais. 


2 — CINÉTICA DO ESCOAMENTO 


À —- Coordenadas Intrinsecas 


No início deste estudo começamos por admitir um 
escoamento sem perdas no interior da roda de uma 
turbomáquina radial com um número de pás infinito. 

Consideremos um referencial xOy, solidário com a 


* Este estudo foi realizado durante o periodo em que o autor prestou serviço no Laboratório de Fisica e Engenharia Nucleares (J E N). 
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roda da máquina, e um referencial XOY, fixo. Utilizando 
um sistema da coordenadas intrínsecas, vamos designar 
por * a 
F = Fr F| (2.1) 


o vector de posição de uma partícula de fluído P; por 
convenção (figura 1) 


r é o módulo do vector r 
r é o versor segundo a direcção de r 
9 é o ângulo de r com o eixo Ox 


q é o ângulo de Ox com o eixo OX 


E imediato que o vector 


(2.2) 


tem também módulo unitário e direcção perpendicular 
A velocidade relativa (no referencial x0y) do ponto 


d ; o | do 
W = na ed 1 el o 
dt dt dt (2.3) 


Derivando novamente obtemos a aceleração rela- 
tiva: 


uma vez que 


do, 


da (2.5) 


A velocidade do transporte (do ponto P suposto 
fixo em xOy) é dada vor 


dt (2.6) 


dé 
dt (2.7) 


obtemos 


U=wr é (2.8) 


A velocidade absoluta (no referencial XOY) define- 
-se por: 


— 


e=u0U-w (2.9) 


ou, por substituição 


* O sinal e simboliza o produto escalar de vectoras 


a dr dg , 
= —— = ti ter + wl 8; 
dr dt (2.10) 
A aceleração de transporte (ou centrifuga), prove- 
niente do movimento de rotação da roda, é: 


== 


dt (2.11) 


o a aceleração de Coriolis a, , perpendicular, em cada 


ponto, à superfície da pá (a, « w = 0)*, com origem na 
acção da pá em movimento sobre o fluído: 


de dr 


= dao = -— 6 
C co ( dt : dt (2.12) 


a 
À aceleração absoluta 
c=w +E' +a, (2.13) 


pode, em face dos resultados anteriores, apresentar a 
forma 


= E: de E d-a 
Cc = de ae o | EoP IREI E o—A 
dt dt À dt” 


dr da E 
+ E— | ——— Je) B; 


(2.14) 


B — Triângulo de Velocidades 


Numa turbomáquina em movimento as trajectórias 
das partículas de fluído podem ser observadas segundo 
dois referenciais diferentes: relativamente ao referencial 
ligado à roda ou ao referencial fixo da máquina; as ve- 
locidades são, também, diferentes em cada um daqueles 
referenciais, designando-se, respectivamente, por velo- 
cidades relativa e absoluta, ligadas vectorialmente, como 
vimos atrás, pela velocidade de transporte. 

Às linhas de corrente, paralelas à velocidade rela- 
tiva, são impostas pela forma das pás da máquina, ge- 
rando-se um campo de velocidades na roda para cada 
velocidade angular de rotação. 

Nem todas as componentes da velocidade são de- 
pendentes da forma das pás: a velocidade de transporte 
só depende da velocidade de rotação e a velocidade 
radial: 

dr 


À dt (2.15) 


para fluidos incompressiveis, é apenas função do cau- 
dal escoado pela máquina. 

Designando por triângulo de velocidades a repre- 
sentação da composição vectorial das velocidades num 
ponto de um escoamento estacionário (figura 2) defi- 


= 


nimos os ângulos 


o] 


a = formado pelos vectores 


mt 


à = formado pelos vectores U e w 
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O ângulo 5 detine a inclinação tangencial das linhas 
dao corrente no referencial solidário com a roda. Conven- 
cionando afectar de um índice r a projecção de um 
vector segundo a direcção radial e de um indice 6 a pro- 
jecção da um vector segundo a direcção tangencial, 
podemos escrever, por aplicação imediata das expres- 
sões das velocidades: 


dr 
E, dt dr 
tous -——— E, 
co de dt 
F(—— + q) 
dt (2.16) 
MEM a 
o we r da ea ed 
De (2.17) tiram-se as relações 
dg cotg f 
“dem r (2.18) 
dºg cosec'5 df cotgf 
do ré dr r (2.19) 


qua serão utilizadas na obtenção de equações particula- 
res da velocidade e aceleração. 

Podem, também, tirar-se da figura 2 algumas rela- 
ções fundamentais, tais como: 


cw, =w seng=c sena (2.20) 


U=G-—-wp=CCcosa-w cos f (2.21) 


Na construção dos triângulos de velocidades usual- 
mento utilizada (figura 2-ll) define-se a posição do vec- 
tor w relativamente ao vector à mediante um ângulo 8, 
relacionado com f por 


B'=]u-—B (2.22) 
podendo então escrever-se (2.21) na forma 
ú=ccosut+w cos pf (2.23) 


C — Linhas de Corrente 


À forma das linhas de corrente no referencial soli- 
dário com a roda, o único tipo de linhas de corrente 
referido nestas considerações, está intimamente relacio- 
nado com a forma e disposição das pás, podendo admi- 
tir-ss uma sobreposição perfeita num escoamento com 
um número de pás infinito, 

Partindo da expressão (2.17) e conhecida a relação 


B = Blr) (2.24) 


para um tipo de pás em estudo, podemos obter, por in- 
tegração: 
ºr cotg B 
À — &4 = | E dr 
A r (2.25) 


ou fazendo 4, = À sendo À constante e OA 2%: 
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*r cotg 8 
Emo Cp e Ke 


Ta (2.26) 


designada por equação geral das linhas de corrente. 

Gerando uma superfície cilíndrica deslocando uma 
recta, paralela ao eixo da roda, sobre uma linha de cor- 
rente (figura 3), obtemos uma superfície que convencio- 
namos designar por superficie de corrente. 

O comprimento elementar ds de uma linha de cor- 
rento obtém-se a partir da aplicação do teorema de Pitá- 
goras a um triângulo infinitesimal (figura 4). Então: 


(ds)' = (r de)" + (dr)º (2.27) 
ou 
ds | dg = 
dr V LC uam 
ou ainda, recorrendo a (2.18) 
ds =-cosec 3 dr (2.29) 
Integrando esta expressão entre r. & r; 
ed 7 
5 = ) : cosec 5 dr (2.30) 


obtemos o perímetro da linha de corrente entre dois 
pontos de posição conhecida. 


D — Equipotenciais 


As linhas perpendiculares em cada ponto às linhas 
de corrente denominam-se equipotenciais, cuja equação 
se obtém a partir da expressão (2.17) uma vez que ao 
ângulo de inclinação tangencial da linha de corrente 
corresponde um ângulo de inclinação tangencial da equi- 
potencial 8 + m/2. Então 


1 dr 
o a o qa (2.31) 
ou 
1 dr 
PRECO o a (2.32) 
integrando obtemos 
“ro tag 
= — | dr + » 
Jr r (2.33) 


em que » é constante e O<vr< 2% 


iImaginemas agora que pretendemos estudar uma 
propriedade de escoamento ao longo de uma normal 
entre duras linhas de corrente. Admitamos que as ex- 
pressões analíticas dessas duas linhas de corrente são 


*r cotg f 
= hr | - == EP 
H R Fi F (2.26-A) 
Pr cotgf 
= aa | aa E 
LE! r (2.26-B) 
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e da normal 


É tg 
p=r,— | dr 
E dm ? (2.33) 


e que a normal encontra a primeira linha de corrente no 
raio r , conhecido, e a segunda no raio r b des- 
conhecido. 


Para r=r, vem 6=t;e 


“r, cotg f 
atá 4 
E Ra dm E (2.34-A) 
iu” tg 8 
ta = "k SE | a 
: Pr r (2.34-B) 


ou, por eliminação de 0, explicitação em Me 


“ra t9B + cotg É 
= hat | Siege care ed 
' es F r (2.35) 
Para r = rp vem j = 04, & 
“r, cotg f 
1,54, + [E a 
F, I (2.36-A) 


CE tg B 
p= o — dit 
b K | Ff, Ff 


(2.36-B) 


donde eluminando agora 9, & substituindo rj pelo 
seu valor obtemos: 


ra t9B+cotg A "ro, to B + cotg 
pe dr= | E. 
r, 


dr = 
. r Rc r 
= Às A (2.37) 
ou ainda 
“Tr, tg B + cotg 8 
| —g == (328) 
fi r 


Depois de feita a integração e resolvendo esta 
equação em ordem a Fp obtemos o seu valor que tinha- 
mos admitido desconhecido, 


3 — COMPONENTES PARALELAS E NORMAIS AS 
LINHAS DE CORRENTE 


A — Transformação de Coordenadas 


Determinaremos agora a forma de decompor um 
vector, representado pelas suas componentes radial e 
tangencial, segundo as direcções da tangente e da nor- 
mal às linhas de corrente. Então, atendendo às conven- 
ções feitas: 


W=W rw, (3.1) 


Introduzindo um verso: z,, paralelo ao eixo da roda 
tal que define com r. e €; um triadro directo (z, X r, = 


= 6,)º obtemos um vector perpendicular a w realizando 


* O simbolo X representa o produto externo de vectores. 


vw produto vectorial Z1 X w ou, atendendo a (3.1): 


aXw= Wo fi + w. 6 (3.2) 


A 
Resolvendo o sistema constituído pelas equações (3.1) 
o (3.2) em ordem a re 8, obtemos: 


wW Wv 


== r a- A E RRSEs 

fã E o A] = EE AM) 
wo w (3.3) 

= Wo a 

di = we + Z X W 
wº o ta (3.4) 


quo nos definem a decomposição dos vectores f, E 4; 


segundo as direcções we lz X w). 


So preferimos utilizar os versores daquelas direc- 
ções 


= | 


ow (3.5) 
nXw 
“mo (3.6) 
podemos escrever 
n=senft-cosfn. (3.7) 
W=cosB tr +seng n (3.8) 


Esta transformação de coordenadas node ser efec- 
tuada por meio de um operador matricial da forma 


sen5 cosf 


-cos 3 sen (3.9) 
Assim, qualquer vector da forma 
P=y GM HS cr 
RS e vs (3.10) 
pode ser transformado à custa do produto matricial 
ne sen cosS 7 
E == ke | — 
-cos 5 sen B vg; 
o vsenf + vcosfã 
—v cos +v sengjf (3.11) 


ou 
v= lv sen + v cos 5) ty + (— v, cos j + 


+v senf)n (3.12) 


E — Velocidade 


Substituindo a relação 
de da dr 


dt = do da. (3.13) 
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nas equações das velocidades obtemos: 


= dr o de 
v=—-——— In+T 84) 

dt dr (3.14) 
Ú=uro (3.15) 
” A do dr 
Cc = + r + w) 6 

dC Cdr dt E (3.16) 
e substituindo (2.15) e (2.18) 
W=w ; (r. + cotg 8 61) (3.17) 
à=oerêo (3.18) 
c= W, n+ tw, cotg A + wr) A (3.19) 


Aplicando a estes vectores o operador (3.9) obte- 
mos a sua decomposição segundo as direcções paralela 
a normal ao escoamento: 


w=cosecfã wt (3.20) 
Ú=uwr (cosfB t+senf ni) (3.21) 
E= (m , cosec 8 + wr cos 8) ti + qr sen 8 n (3.22) 
C — Aceleração 
Derivando (3.13) em ordem ao tempo 
d'a da dr do dr 
SR! ed Ria dr 
dt dr dt dr dt” (3.23) 


o substituindo (3.13) e (3.23) na equação (2.4) da ace- 
leração relativa obtemos: 


o d'r dg dr 
“= — ei ff À + d=ero lh 
= dt” ar! a | 
dr de dr dg 
+ 2 ( e = Es cpa + r | = Jº = 
dt dr dt dr 
da Er 
“= p a i À 
dr dr (3.24) 
ou atendendo a (2.18) e (2.19) 
| dºr cotg' & EA 
w' = —— - woln+H|l w. + 
| dt r ! r 
cotg £ dB dr 
(| mi = re — ) + cotg £ 9 
r p dr 9 é dt | 
(3.25) 


As equações das acelerações relativa e de trans- 
porte podem também escrever-se 


É or rr (3.26) 


Es E Zow, | — cotg 8 nr + 64) (3.27) 


Cl 


Compondo (3.25), (3.26) e (3.27) através de (2.13) 
obtemos: 
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tr + 


dt r 


= | dºr tcotg Bw, + or) 
| == —=6 s 


w, | dp 
ri ——— tw cotg 8 + 20r) — wº. cosec” 8 F + 
r Fr 


dr - 
de |” 


Aplicando novamente o operador (3.9) aos vecto- 
res (3.25), (3.26) e (3.27) obtemos as componentes 
das acelerações no referencial definido pelos vectores 


+ cotg 5 
q (3.28) 


ti B ny 
” | dºr dB o 
w' = cosec À | — wº cot Es: É 
Í , dt? g B dr str 
cotg É dj e 
de E a nm 
' r dr (3.29) 
u=owr(-senfB t+cosB n) (3.30) 
aç = 2 ww, cosec f n (3.31) 


a por adição das equações (3.29), (3.30) e (3.31) 


E E dºr . dB 
c' = cosec f | ( — — W, cotg 8 pés 
dt” dr 
| e — cotg B . dB 
= E sa) AT =, FF efe 


co" F gas 
E 2 P 
+- 5 sen 28 +2ww,) n | (3.32) 


Os vectores velocidade aparecem-nos dependentes 
da velocidade radial w , imposta pelo caudal, ângulo 
característico da pá e velocidade angular e os vectores 
das acelerações dependem, além dos valores atrás indi- 


d'r 
cados, da aceleração radial dê e da derivada e - 
x Fr 


q — ACÇÃO SOBRE AS PAS 
A — Efeitos da Aceleração 


As diferentes parcelas que compõem a aceleração 
absoluta têm efeitos diversos na distribuição das pres- 
sões no interior da roda (figura 5). 

A aceleração relativa pode ser estudada admitindo 
a roda parada com escoamento; o fluído, ao descrever 
o canal curvo entre pás muda de direcção, gerando um 
gradiente do pressões, dirigido para a parede côncava 
do canal. Os efeitos da aceleração de transporte podem 
ser postos em evidência admitindo a rotação da roda 
sem escoamento: as pás sofrem uma pressão radial, di- 
rigida da dentro para fora, e devida à aceleração cen- 
trifuga do movimento. Finalmente, a aceleração de Co- 
riodis, proveniente da acção directa das pás sobre o 
fluído, impele-o contra a pá seguinte no sentido do mo- 
vimento. 


ol 


A distribuição de pressões resultante da actuação 
simultânea das três componentes pode ser identificada 
com a distribuição devida à aceleração absoluta. O es- 
tudo feito conduziu à decomposição deste vector acele- 
ração absoluta segundo as direcções paralela e perpen- 
dicular ao escoamento, admitido como unidimensional; 
a componente paralela influencia a variação de pressão 
ao longo da roda e a velocidade de escoamento, en- 
quanto que a componente normal gera uma diferença 
de pressão entre as faces das pás. Este último fenómeno 
será estudado a seguir. 

Bb, análise do efeito de um escoamento num canal 
poda fazer-se através do teorema de quantidade de mo- 
vimento, enunciando que uma partícula de fluido de 
massa dm, sofrendo uma variação instantânea de velo- 
cidade de, processa sobre as paredes do canal uma 
acção dF aplicada durante o intervalo de tempo dt tal 


que 
dF dt = dm de (4.1) 
ou ainda 
" dm 
ida ato fe (4.2) 


dm 


Como a fracção reoresenta o caudal mássico G 


do escoamento, obtém-se, por integração 
E = G (cs — c,) (4.3) 


Os valores c; e c são tomados entre as duas Secções 
quo delimitam a zona do escoamento cuja acção pre- 
tendo ser estudada. Se aplicarmos o teorema de quanti- 
dado de movimento entre as secções de entrada e saída 
da roda de uma turbomáquina obtemos a acção total do 
escoamento sobre as pás, 


B -— Pressão Normal 


Quando se pretende calcular a diferença de pressão 
entro as faces de duas pás consecutivas é necessário 
deixar de admitir um escoamento unidimensional pois 
que o fluído só encontra distribuições de pressão análo- 
gas se percorrer linhas de corrente com posições simi- 
lares relativamente às pás. 


Nestas condições temos de admitir um escoamento 
bidimensional em que todas as particulas de fluido per- 
correm trajectárias semelhantes, 


Consideremos duas pás consecutivas suficiente- 
mente próximas e afastadas angularmente de ds, e 
dois planos perpendiculares à superficie das pás e afas- 
tadas de uma distância elementar ds. Se as pás tiverem 
uma largura |, o elemento de volume dV (figura 6) li- 
mitado pelas áreas dS sobre a superfície da roda 
será: 


dV — | dA (4.3) 


ou, sendo | o ângulo que a pá forma com a direcção 
tangencial 


2 


dA = r dr de (4.4) 


As áreas dS | e dS, são iguais a 


dr 
da ="— | 
nº senf (4.5) 
dS, = r de sen B | (4.6) 


Sendo p a densidada do fluido em escoamento, a 
massa dm do elemento de volume dV será: 


dm = plr dr de (4.7) 


o admitindo dp a A diferença d epressão entre as su- 
perfícios dS, e dp, a diferença de pressão entre as 
superfícies dS, , as forças de pressão aplicadas ao vo- 
lume elementar dp, ds, e do, ds, deverão equili- 
brar a força de enércia. Então 


dp,dS,n + dp, dS, t+dm c'=0 (4.8) 
ou ainda, igualando as componentes homólogas 
dp, dS,m= —plr dr do e, (4.9) 


dp, dS, t= —plr dr de c, (4.10) 


ou atendendo a (4.5) e (4.6): 
dp,= —prf senf c, dy (4.11) 


dpt = — p cosecf c', dr (4.12) 


à diferencial do, representa a variação elementar 
de pressão ao longo da linha de corrente, enquanto dp. 
é a diferença de pressão infinitesimal entre duas pás 
consecutivas. 

A integração das equações (4.11) e (4.12) conduz 


Ap,=—P | r sen f Cc, do 


) (4.13) 
APp,=— | cosec /) c, dr 


(4.14) 


À segunda integração é imediata dentro da hipótese 
de unidimensionalidade admitida e será efectuada 
adiante. 

A primeira integração exige o conhecimento do 
número de pás. 

Um caso existe em que a integração é imediata: 
quando o percurso segundo o qual ela é efectuada é um 
arco de circunferência. Então os valores de à er são 
constantes e, sendo z o Número de pás da roda, o afas- 
tamento angular 6, entre duas pás consecutivas será 


id z (4.15) 


Integrando então (4.11) ao longo de um arco de 
circunferência obtemos 


à p= — prseng co, 


(4.16) 
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C — Pressão Tangencial 


Para determinar a pressão tangencial ao longo de 
uma linha de corrente substituimos em (4.14) o valor 


dr — dB 
G.= coseo | ( rodo w'. cotg 8 = 
= jd sen" 8) ) (4.17) 
obtido por (3.32) ou: 
- dºr 
dp, = — p cosec o( re 
dB 
— w cotg 8 — — o, sen" f ) dr (4.18) 
ou, integrando 
* Fã dºr 
4 Dj = — | (cosec” f a - 
; dj E 
= w', cotgB cosec' : — Cego) dy (4.19) 


a desdobrando o integral numa soma de duas parcelas 


Ps dºr dr 
A pç=— | N (cosec' 8 o dt 
dj “Fa 
Di dd PRA AT 


A primeira parcela pode escrever-se 


"tl a*r dr 
cosec ——— potg É cosec” 5 
| É | P dt” dt 9 É f 
da |] dr Pts dr dr 
peca doce lies | cosec 8 e — aim 
dt dt J t dt dt 
dr dB | 1 
— JH cotgfcosecp-——— Idt= 
dt ada Pede | 2 
dr Tt 1 
too a —— == = | W— Wº, 
dt to 2 z 
(4.21) 


A segunda parcela é de cálculo mais simples 


Finalmente, por substituição em (4.20); 


| (w'— u5) e (w | = a? | (4.23) 


A 2 


5 — ALTURA TEÓRICA DE ELEVAÇÃO 
A — Equação de Bernoulli (Fluidos Incompressíveis) 


A altura da elevação de uma turbomáquina defi- 
ne-so como sendo a energia transmitida à unidade de 
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peso do fluído em escoamento posta em jogo sob duas 
formas diferentes: energia cinética e energia de pressão, 
correspondendo a cada uma delas uma parcela na equa- 
ção da altura teórica de elevação. Designando por: 


p— a pressão num ponto do escoamento 
g—a aceleração de gravidade 


Yy— o peso específico do fluido suposto incompres- 
sível 


e distinguindo pelos índices 1 e 2 dois pontos do es- 
coamento respectivamente à entrada e à saída da roda, 
a altura de elevação H adquirida pelo fluído no interior 
da roda será: 


29 a no (5.1) 


já conhecida como equação de Bernoulli aplicada a um 
escoamento sem perdas. 

A energia cinética adquirida pela unidade de peso 
do fluído é dada pelo termo 


= 
Cs C4 


29 


O valor do acréscimo de pressão verificado no 
fluido durante o seu escoamento no interior da roda 
poda determinar-se a partir da composição de dois mo- 
vimentos estacionários e adicionando algebricamente o 
aumento de pressão devido cada um deles isoladamente, 
aplicando o princípio de sobreposição dos estados de 
equilibrio. 

Designando por p:— p; o acréscimo de pressão 
devido ao movimento de rotação da roda sem escoa- 
mento e por p" —p” o acréscimo de pressão criado 
pelo escoamento do fluído admitindo a roda parada, 
temos: 


P, Esc P, = (p'. A p'.) E (p”. as p“,) (5.3) 


Pela definição dos dois efeitos é fácil de prever 
que a primeira parcela de (5.3) dependerá da veloci- 
dade de transporte enquanto que a segunda será função 
da velocidade relativa do escoamento ou, de outra forma, 
a primeira depende da velocidade de rotação da roda 
a a segunda está relacionada com a curvatura das pás. 
Faremos a seguir uma análise da cada uma destas par- 
celas. 


B — Rotação sem Escoamento 


Admitamos que uma roda de uma turbomáquina ra- 
dial roda sem escoamento e consideremos no seu inte- 
rior duas superfícies cilíndricas coaxiais distanciadas 
de dr. Designemos por dS, e dS. duas áreas elementares 
definidas sobre as superfícies cilíndricas por um diedro 
de aresta coincidente com o eixo e amplitude do, e re- 
presentemos por dV o volume infinitesimal limitado por 
aquelas duas áreas (figura 7). Então, para uma roda de 
largura |; 


ER: 


dS,=I rode (5.4) 
dS.= | (r+ dr) do (5.5) 
dV =| Fr do dr (5.6) 


A força centrífuga F. a que esta sujeito o elemento 
dV deverá ser: 


Fe=p dv or (5.7) 


em que p é a massa específica do fluído. 

A força devida à pressão infinitesimal dp' deve 
equilibrar a força centrífuga originada pelo movimento 
de rotação: 


p dV w r= dp' ds (5.8) 


ou, por substituição de (5.4), (5.5) e (5.6) e despre- 
zando infinitésimos de ordem superior 


wu 
dtpri=a» Ff dt= = dr 
p pu É ; (5.9) 
Como 


Y=p9 (9.10) 


integrando (5.9) entre os pontos 1 e 2 obtemos: 


Ff 
* Em é ap - | tr = r*) 
Bb P| 29 u ? Ni) 
ou ainda, por (2.8): 
Ê r + ( a : ) 
.— E O As = a 
di 29 (5.12) 


U— Vóártices 


Quando um escoamento percorre um canal curvo 
a variação direccional da quantidade de movimento do 
fluído provoca uma acção sobre a conduta; este efeito 
é a consequência do gradiente interno de pressões 
criado pela presença da componente centrífuga da ace- 
leração. Supondo um escoamento sem perdas, o princi- 
pio da conservação de energia permite afirmar que, não 
havendo trocas energéticas do fluído com o meio, a velo- 
cidade das partículas sujeitas a maior pressão terá que 
diminuir a fim de que a queda de energia cinética Es 
conpense o aumento de energia de nressão E, sofrida 
pela unidade de peso do fluído, ou: 


E e E, = Constante (5.13) 


Admitiremos que se o escoamento se der numa 
conduta de curvatura circular (figura 8) as linhas de 
corrente são arcos de circunferância*, Designando por 
u, a velocidade do fluído à entrada da curva e por u; a 
velocidade sobre um arco da linha de corrente à distân- 
cia r do centro, a condição (5.9) permite escrever 


dp 1 Tm 
; E = dr 
t q r (5.14) 


+ Esta hipótese não é rigorosa para um fluido real 


4 


& diferenciação da expressão (5.2) da energia 
cinética conduz a 


u du 
de —>—— 
p g (5.15) 
Diferenciando (5.13) 
dE, + dE, —O (5.16) 


a substituindo (5,14) e (5,15) 


e (5.17) 


u, en (5.18) 
podendo inferir-sa que 
u r=> Constante (5.19) 


Um escoamento com estas características toma O 
nome dos vórtice e, quando verificado numa turbomá- 
quina, traduz a ausência de energia comunicada pelas 
pás ao fluido. 

Para escoamenitos de curvatura não circular é 
legitima a determinacão do raio da circunfrerência oscu- 
latriz correspondente a um pequeno troço da linha de 
corrente a fim de determinar, por aproximação, a distri- 
buição de pressões local. 


D — Escoamento sem Rotação 


Vamos sulpôr que o escoamento se dá encon- 
trando-se a roda da turbomáquina parada. Uma vez que 
não há energia fornecida ou cedida ao escoamento a 
aquação dae Bernoulli toma a forma 


Par pos Wºs = Wr, 
Y 29 (5.20) 
ou ainda 
Y 
Do— po = = (W', — Wº,) 
2g (5.21) 


E — Rotação com Escoamento 


O valor de p; — p. é obtido a partir das relações 
(5.3) (5.12) e (5.21): 
? E | 
= p= sus) AE ua] | 
29 : (5.22) 


a substituindo esta expressão em (5.1) obtemos 


H = = Peg = 60) + tuto — Ui j— imo — 2] 


[5.23] 
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A análise do triângulo de velocidades (figura 2) per- 
mita escrever 


co+u-w = (e dede te, + mm = qu WO 
(5.24) 

ou, simplificando 

Cro -m=2to wa] é: (5.26) 

ou ainda 

e +u-w=2uc (5.26) 


f 


Substituindo esta expressão na equação (5.23) 
obtemos ainda 


H = (ué Ch, — Ui Ch) (5.27) 


expressão simplificada da altura teórica de elevação de 
uma turbomáquina radial. 
O produto u cg pode ainda ser desenvolvido: 


ucy=ulu—w)=ulu-— w, cotg B) (5.28) 
donde, por substituição em (5.27): 
1 
H = fu”; — uS) + (us W cotg B: — Ui W 
g : 
cotg &, 
Bd | (5.29) 


6 — BINÁRIO E POTÊNCIA 
A — Binário Resistente 


Idealizemos um escoamento entre duas pás muito 
próximas de uma turbobomba radial em que 1 e 2 são, 
respectivamente, dois pontos à entrada e à saida do 
tubo da corrente definido pelas pás. Aplicando o teorema 
do momento angular em cada um dos pontos referidos, 
obtemos: 


(6.1) 


dt i a (6.2) 


em que m é a massa do fluido contido no tubo de 
corrente, No ponto 1 o fluido começa a ser arrastado 
pela roda originando um acréscimo de binário resistente 


M,. Por convenção 
AM >0 (6.3) 


No ponto 2 o fluido deixa de sofrer a acção da 
roda dando-se uma diminuição de binário resistente M., 
implicando 


O binário resistente M será então 


NOTA: A equação (5.22) já tinha obtido por outra via em (4.31) 
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M=— AM — à Ms (6.5) 


As equações (6.1) e (6.2) foram obtidas para pon- 

tos definidos do escoamento, verificando-se aí que 
dr 

dt E (6.6) 


admitindo o escoamento estacionário, pode também 
escrever-se 


de 
—.... = () 
dt (6.7) 
concluindo-se então que 
E ER! =x ms 
AM=— r Xe 
dt (6.5) 
Decompondo o vector 
Cc=€c +c (6.9) 


obtemos, por substituição 


AM=Grco a (6.10) 


em que 
z, é o versor paralelo ao eixo da máquina 
G é o caudal mássico 
Retomando a equação (6.5) podemos escrever, 


finalmente: 


M=6 tt ur, 0,) (6.11) 


B -—. Bombas e Turbinas 


Consoante uma máquina funciona fornecendo ou 
recebendo energia assim o binário é motor ou resistente, 
tomando respectivamente o sinal positivo ou negativo, 
por convenção. No caso particular das turbomáquinas 
temos: 

M > O para bombas 
M <!O para turbinas 


ou, recorrendo à equação (6.11) 


fz Cu > [4 Chy para as turbinas 
rs Ch. < Mm Coypara as bombas 


É, assim, fácil de concluir que a solução constru- 
tiva mais provável indica 
rs > rf; para as bombas 


rs < r, para as turbinas 


Por esta razão nas turbomáquinas radiais, salvo 
raras excepções, as linhas de corrente das bombas são 
descritas do centro para a periferia da roda enquanto que 
nas turbinas o sentido do escoamento é inverso, 


C — Potência 


Comparando as expressões (5.27) e (6.11) con- 
cluímos que, em valor absoluto, a potência P absorvida 
ou cedida pela máquina é: 


P=Muwu=goGMH (6.12) 


Consoante o tipo da turbomáquina assim 


P > O para as bombas 


P << O para as turbinas 


7 — ESCOAMENTOS RADIAIS INCOMPRESSÍVEIS 
à — Velocidada e Aceleração 


Sendo O o caudal de fluido, suposto incompres- 
sivel, escoado por uma turbomáquina radial, e repre- 
sentando por | a largura útil da roda, temos que 


dr 
dt (7.1) 


Q=2g7ril 


Admitindo um caudal invariável, podemos introduzir 
a constanta 


Q 
age > 
2 q | (7.2) 
e portanto 
dr K 
dr or (7.3) 


Derivando (7.3) em ordem ao tempo obtemos 
dr uy" 
dr r (7.4) 


Substituindo (7,4) nas equações das acelerações 
(3.25) à (3.28) obtemos: 


, mr o dB 
= — ses bio is = Pa) 
r dr (7.5) 
| 1 dh q 
e eCrnTirwW----=tuw.) & 
r dr F 


(7.6) 


A decomposição destes vectores segundo as direc- 
ções tangente e normal às linhas de corrente pode fazer- 
-se utilizando o operador (3.9) 


w E 
w' = - — | (seng +r cosp-——)t, + (r sen A. 
r dr 
dB o 
E cos gn, 
dr Ê | (7.7) 
a sen d | dB 
c' c +reotg A w 2 
| gw a UW,. 
dB see 
cotg 8) ty, + (r w ” 2 uw, — € cotg B) m | 
r 


(7.8) 


O ânguloa, definido pelos vectores c e u, pode 
calcular-se recorrendo a (2.16] 


u 


cotga — + cota À 
e” w (7.9) 


8 — Pressão sobre as Pás 


Admitindo um escoamento unidensional com um 
número de pás suficientemente grande, a diferença de 
pressão entra as duas faces de uma pá num ponto é 
igual à diferença de pressão, no mesmo canal, entre 
duas pás consecutivas em pontos à mesma distância 
do eixo. 

Substituindo em (4.16) a componente normal da 
aceleração absoluta obtida em (7.8) 


2 —- dB 
Ba = ” e semp rw E 2 u Wi 
— jp? cotg £) (7.10) 


Efectuando a integração (4.14) segundo o mesmo 
caminho (r constante) temos que 


Ap, =0 (4.4] 


ou seja a pressão tangencial mantém-se constante, e as 
pressões num ponto da parede de uma pá são 


prep, + áàp, (7.12) 


C — Altura Teórica de Elevação 


O desenvolvimento da equação (5.29) conduz a 


a a 
H=——— | us— ui +oK (cotg 8 — cotg 81) 


(7.13) 


No caso mais favorável de entrada (ausência de 


pré-rotação) o ângulo dos vectores ceu deverão ser 
a => 7/2 donde 


LU 


cotg 8, — — ; 

Cc, (7.14) 

Como nestas condições a velocidade absoluta só 

tem componente radial é c; =W e por substituição 
em (7.13): 


1 


— tus + us w. cotg 8.) 
g ty SESRE (7.15) 


3 E 
|| 


Introduzindo a constante (7.1) na equação anterior 
podemos escrever: 


L”3 uu 
=—— cotg 3. OQ 
g é | 9 É (7.16) 


a ii 


equação que se designa por característica teórica de uma 
turbomáquina com entrada radial (figura 10). 
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8 — ESCOAMENTO A ENERGIA CONSTANTE 
A — Características 


Designamos por +. scoamento neutro aquele em que 
a energia total do fluído se mantém invariável durante 
todo o percurso na roda. Então, a energia posta em 
jogo pela máquina é nula, assim como a potência e o 
binário. 

O anulamento da expressão (5.27) obriga a que 


UC; =u; Cy, (8.1) 
Ou, numa outra forma 
FC CcOSasz tc cosa; (8.2) 


As relações entre os elementos geométricos do triân- 
gulo de velocidades (figura 2) permitem escrever 


E 
E (8.3) 


Cc cos a = V E — 


e introduzindo esta relação em (8.2) 


r Pta a . = f e, ue Ee 
Ve- SM z ga 
ou simplificando 
Cr=€ (8.5) 


Comparando as equações (8.2) e (8.5) verificamos 
que «e = a, OU seja, o ângulo « mantém-se constante ao 
longo da linha de corrente. Admitamos então 


E — je (8.6) 


Quando w= O então c>=0,24=0 820 movimento 
degenera em vártice, não havendo escoamento radial. 
Quando «= 7/2 vem c= C, - podendo escrever-se 


(8.7) 


a o fluído efectua um escoamento puramente radial ou 
fonta em que o movimento da pá não intefere. 

Os vortices e as fontes serão analizados como 
casos particulares dos escoamentos a energia constante, 


C=uW+w 


B — Velocidade e Aceleração 


Uma vez que o ângulo « é constante ao longo do 
escoamento a relação (7.9) permite-nos escrever 


L 
cotg 8 = cotga — ———— 
a W, (8.8) 
ou por (7.3) 
dd 
cotg 3 — cotga | — K r (8.9] 


Derivando esta equação em ordem a r a atendendo 
ainda a (7.3) obtemos 


df ui 


dr wW r (8.10) 
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Introduzindo as relações (8.8) e (8.12) em (7.5) 
a (7.6) podemos escrever: 


1 = = 
SM FZ Ud) 
Fr 


(8.11) 


wo. 
c' == - Coser a 
r H 


(8.12) 


A aceleração absoluta apresenta apenas compo- 
nente radial, condição esta que permite confirmar o anu- 
lamento do binário e potência em jogo na roda, pois que 
a resultante das forças aplicadas sobre um elemento do 


fluído deverá ter a mesma direcção que c. 
Particularizando temos, para o caso do vórtice 


w' = 0 (8.13) 
o =0 (8.14) 
a para o caso da fonte, atendendo a (8.7) 
" TO = 
Ww = — aci (Cc nt2uw 6) (8.15) 
Gero 
E O (8.16) 


C — Linhas de Corrente e Pressão sobre as Pás 


As equações (2.17) e (8.9) permitem escraver 


dá uw 


p= = .C01 = I 
dr ME (8.17) 


ou, separando variáveis e integrando, admitindo que 
para r=r, é 8=0 


r qo 
e=cotge In -— — E = E] 
A r, 2K (8.18) 
ou ainda 
r 1 u Us 
8 = cotg An -— em — — ) 
E; 2 W W 1 
(8.19) 


equação que caracteriza as linhas de corrente. 

Para um vortice, como a = O é necessário que 
r=r, para é ter um valor finito, concluindo-sa então 
que para este caso a equação das linhas de corrente 


será tal que 
t = constante (8.20) 


Para uma fonte é E = q/2 e então 


1 u U; 
2 = E e (8.21) 


F 


q = — 


À pressão sobre as pás da máquina pode determi- 
nar-se substituindo a relação (8.10) em (7.10): 


2% 
» Sen B cos B c 
o O PRE RO (8.22) 


AP, =— 


7 


9 — ESCOAMENTOS COM AUMENTO DE ENERGIA 
A — Velocidade Relativa Constante 


Admitindo que a velocidade relativa ao longo do 
escoamento se mantém constante 


PEÇA (9.1) 
vem 


W, K 
sen 8 =- = — 


wi War (9.2) 


Derivando esta expressão em ordem a r podemos 
escrever 


dp tgB 


dr r (9.3) 
a substituindo em (7.5) e (7.6) 


w = —— (nr — tg8 8:) 


(9.4) 


1 


(9.5) 


| e n-(Zu W + wo taof) 9; 
r LA 
Decompondo w' segundo as direcções tangente e 
normal às linhas de corrente recorrendo à equação (7.7) 
verificamos que 


wr 
F, 


|| 


w' “secj n 


(9.6) 


concluindo-se que a aceleração relativa é normal, em 
cada ponto, à linha de corrente. 
A equação (2.18) pode escrever-se 
de cotg 3 dr 
dB r dA (9.7) 


Substituindo nesta equação o valor obtido em (9.3) vem 


dg 
ET a id (9.8) 
obtendo-se, por integração, a equação das linhas de 
corrente 
6 = (cotg 8 — cotg 8:) + (8 — Bi) (9.9) 


O perimetro de um segmento de uma linha de cor- 
rente, limitada por dois pontos à distância r e r; do 
centro da roda, obtém-se substituindo em (2.30) o valor 
(9.2) e integrando 


Wi 
tr mas: E) 
2K (9.10) 


a ES 
À pressão sobre as pás determina-se introduzindo 
em (7.10) o valor obtido em (9.3) 


27 
Ap, =— Eis: senf (ce cotgB+2uw +w tgB) 


(9.11) 


B — Inclinação Radial Constante 


Partindo da equação (2.17) e admitindo constante 
a inclinação radial, sendo 


B=B, (9.12) 
obtemos imediatamente 


dg cotg f, 


E 


| dr r (9.13) 


Integrando esta equação entre um ponto de raio 1, e um 
ponto genérico de raio r, admitindo que para r=r, é 
9 = O, obtemos a equação das linhas de corrente: 

r 


Fi (9,14) 


9 = cotg 8. Im 


O perímetro de uma linha de corrente entre r e r; 
pode ser determinado por (2.30) que, depois de inte- 
grada conduz a 


s = cosec 5, lr—r1) (9.15) 


Um caso particular deste é aquele em que 5 = 7/2, 
caso das pás radiais planas. Então a equação das linhas 
de corrente transforma-se em 


g = Constante (9.16) 
e O seu comprimento 

= p— E (9.17) 

Pelo facto de ser 


dj 
= O 
dr (9.18) 


as equações (7.5) e (7.6) podem apresentar a forma 


r (9.19) 


(9.20) 


Para as pás radiais planas estas equações transfor- 
mam-sa em: 


= W, 
W = — — E 
r (9.21) 
= Wi - u mi 
Cc =- n—2 84) 
r (9.22) 


F 


O cálculo da pressão sobre as pás faz-se à custa 
da equação (7.12) que pode escrever-se, atendendo a 
(7.10): 


Mer 
Ap, =— === senp,lc' cotg5 +2uw,) 


(9.23) 
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8 — ESCOAMENTO A ENERGIA CONSTANTE 
A — Caracteristicas 


Designamos por « scoamento neutro aquele em que 
a energia total do fluído se mantém invariável durante 
todo o percurso na roda, Então, a energia posta em 
jogo pela máquina é nula, assim como a potência e o 
binário. 

O anulamento da expressão (5.27) obriga a que 


U Cy = Cy (8.1) 
Ou, numa outra forma 
FC COSa=fl CG COSa, (8.2) 


As relações entre os elementos geométricos do triân- 
gulo de velocidades (figura 2) permitem escraver 


Ra 
r (8.3) 


C COS «= Ve- 


e introduzindo esta relação em (8.2) 


Da Kº / Ke 
É = / 
i Er pasa — Fi / c q e 
V r y r (8.4) 
ou simplificando 
Ert=€e ht (8.5) 


Comparando as equações (8.2) e (8.5) verificamos 
qua a = a, OU seja, O ângulo « mantém-se constante ao 
longo da linha de corrente. Admitamos então 


(8.6) 


(É = tp 
Quando w= O entãoc=u,a=0eomovimento 
degenera em vártice, não havendo escoamento radial. 
Quando « = 7/2 vem c = b. | podendo escrever-se 
(8.7) 


a o fluído efectua um escoamento puramente radial ou 
fonte em que o movimento da pá não intefere. 

Os vortices e as fontes serão analizados como 
casos particulares dos escoamentos a energia constante. 


= y+w 


B — Velocidade e Aceleração 


Uma vez que o ângulo « é constante ao longo do 
escoamento a relação (7.9) permite-nos escrever 


u 
cotg 8 = cotga — —— | 
k : (8.8) 
ou por (7.3) 
id 
cotg ) = cotga  — K pé 18.9] 


Derivando esta equação em ordem a r e atendendo 
ainda a (7.3) obtemos 


df u 
(8.10) 


dr wW r 
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introduzindo as relações (8.8) e (8.12) em (7.5) 
a (7.6) podemos escrever: 


wW=--—— (WEt+2uw.6) 
' r (8.11) 


w? 


sms 


cosec” En 


T = -— 
(8.12) 

A aceleração absoluta apresenta apenas compo- 
nente radial, condição esta que permite confirmar o anu- 
lamento do binário e potência em jogo na roda, pois que 
a resultante das forças aplicadas sobre um elemento do 


fluído deverá ter a mesma direcção que c”. 
Particularizando temos, para o caso do vártice 


w' = O (8.13) 
e =" (8.14) 
o para o caso da fonte, atendendo a (8.7) 
Ee É ud a 
“= — — lc n+t2uw 6) (8.15) 
Musa 
r (8.16) 


C — Linhas de Corrente e Pressão sobre as Pás 


As equações (2.17) e (8.9) permitem ascrever 


da FR] 
= cota. —-———— 
dr a 8.17) 


ou, separando variáveis e integrando, admitindo que 
para r=n é0=0 


r o 
o = cotga, In E FR seu 2K E — 4) (8.18) 
ou ainda 
r 1 u aj 
a a Og We Wa ) 
(8.19) 


equação que caracteriza as linhas de corrente. 

Para um vortice, como do E O é necessário que 
r=r, para 9 ter um valor finito, concluindo-se então 
que para este caso a equação das linhas de corrente 


será tal que 
r = constante (8.20) 


Para uma fonte é Rs = w/2 e então 


1 u Ui ) 
2 2 W. W 1 (8.21) 


E — — 


A pressão sobre as pás da máquina pode determi- 
nar-se substituindo a relação (8.10) em (7.10): 


2a” 
sen B cos f c 
7 E i Ê (B.22]) 


a ra 


7 


9 — ESCOAMENTOS COM AUMENTO DE ENERGIA 
A — Velocidade Relativa Constante 


Admitindo qua a velocidade relativa ao longo do 
escoamento se mantém constante 


W =w [9.1] 


Wa (9.2) 


wW K 
sen p = . 


Derivando esta expressão em ordem a r podemos 
escrever 


dB tgp 


dr r (9.3) 
a substituindo em (7.5) e (7.6) 


do 


tn — tgf 61) 


é | 
[ 


(9.4) 


(9.5) 


se 1 Ta ' 
c' = -— [e n—t2u We wo tgB) 94 
r KH 

Decompondo w' segundo as direcções tangente e 
normal às linhas de corrente recorrendo à equação (7.7) 
verificamos que 


wW = 


sec fi n; (9.6) 


concluindo-se que a aceleração relativa é normal, em 
cada ponto, à linha de corrente. 


A equação (2.18) pode escrever-se 


de cotg 5 dr 


dB ro dh (9.7) 


Substituindo nesta equação o valor obtido em (9.3) vem 


de 
= — cotg'5 
dB do (9.8) 
obtendo-se, por integração, a equação das linhas de 
corrente 


o = (cotg 8 — cotg 8;) + (B — Bi) (9.9) 


O perímetro de um segmento de uma linha de cor- 
rente, limitada por dois pontos à distância r e r; do 
centro da roda, obtém-se substituindo em (2.30) o valor 
(9.2) e integrando 


K , 
, (ré - [1] 
2K (9.10) 


SG FE 
A pressão sobre as pás determina-se introduzindo 
em (7.10) o valor obtido em (9.3) 


2% 
AP, =— Es senf (c' cotg8 +2uw tw tgf) 


19.11) 


8 — Inclinação Radial Constante 


Partindo da equação (2.17) e admitindo constante 
a inclinação radial, sendo 


BP =B. (9.12) 
obtemos imediatamente 


da cotg £ 


e 


dr r (9.13) 


Integrando esta equação entre um ponto de raio r, e um 
ponto genérico de raio r, admitindo que para r=r, é 
g = O, obtemos a equação das linhas de corrente: 


F 


F; (9.14) 


o = cotg B, Im 


O perimetro da uma linha de corrente entre r e r, 
poda ser determinado por (2.30) que, depois de inte- 
grada conduz a 


s=cosecf, lr—ri) (9.15) 


Um caso particular deste é aquele em que 4 = 7/2, 
caso das pás radiais planas. Então a equação das linhas 
de corrente transforma-se em 


4 = Constante (9.16) 
e o seu comprimento 

= EE (9.17] 

Pelo facto de ser 


df 


= () 
dr (9.18) 


as equações (7.5) e (7.6) podem apresentar a forma 


mr 


r (9.19) 


wo = —— 


e 1 e se 
G=-=--—[n-Zumw. Bl 

r | r (9.20) 
Para as pás radiais planas estas equações transfor- 


mam-se em: 


= Wo 
wu” a e Conga 
r (9.21) 
” e mo 
E = =—=——=is ln — 2 “8i) 
r w (9.22) 


O cálculo da pressão sobre as pás faz-se à custa 
da equação (7.12) que pode escrever-se, atendendo a 
(7.10): 


Ap, = RS: sen'B le cotg ii, +2uw,) 


(9.23) 
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C — Pã Plana 


Admitamos agora uma pá plana, figura 12, Sejam 
P a P' dois pontos muito próximos sobre a pá, para os 
quais o ângulo de inclinação tangencial da pá toma, res- 
pectivamente, os valores 5 e 8 + df, definindo com o 
centro O da roda um ângulo 


| POP = dg 
Então, pela figura, obtemos 

A+ da = 8 + de (9.24) 
donde 

dg = da (9.25) 
ou integrando e admitindo 4,= O 

p=8— A; (9.26) 


Esta relação, substituída em (2.17) permite es- 
crever 


dB cotg 8 
dr Cr (9.27) 
ou integrando entre r e r; 
a cos à 
In - eia 
F; cos f; (9.28) 
ou, finalmente 
r cosB =" cosgf, (9.29) 


Relacionando o ângulo 5 com 4 substituindo (9,29) 
em (5.26) obtemos a equação das linhas de corrente 


r; 


= are cos (——— cos 8) — & 
r 


(9.30) 
Substituindo (9,27) em (2.30) e integrando ob- 

temos 

” 


sec B df 


5s= nr cos j. | 
Sh (9:31) 


Ou 


= n cos 8, (tgB — tg8.) (9.32) 


equação que nos permite obter o perimetro de um seg- 
mento de uma linha de corrente. 

Substituindo (9.31) em (7.5) e (7.6) obtemos as 
equações das acelerações: 


= wu” — — 
wW = — -——= (1; + cota 8:) = 
a 9 É (9.33) 
= 1 ED = 
E = e co n+iw cotgB-2Zu w.) 6:] | 
Fr 
(9.34) 


A aceleração relativa é, neste caso, paralela à velocidade 
relativa. 
Substituindo (9.31) em (7.12) obtemos 


2 7 u” 
Pp, E 


sen 258 + 2u w.) 
i n (9.35) 
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expressão que nos dá a diferença de pressão entre dois 
pontos da mesma equipotencial sobre duas pás conse- 
cutivas próximas. 


D — Exemplo Numérico 


À fim de comparar as caracteristicas mecânicas de 
diferentes tipos de escoamento vamos admitir que pre- 
tendemos, com uma turbobomba, vencer uma altura de 
elevação de 20m com um caudal de água de 25 |/s. 
As dimensões, admitidas para o impulsor por compara- 
ção com outras do mesmo tipo, [Stepanoff, Centrifugal 
and Axial Flow Pumps, Wiley 1957] são as seguintes: 


Escoamento radial 

Entrada normal, sem choque 

Raio exterior do impulso — 12,5 em 

Raio interior do impulso — 5 cem 

Largura dos canais — 2 em 

O número de pás admitido, a fim de poder manter 
as hipóteses de unidimensionalidade, foi bastante supe- 
rior âquilo que é habitual em máquinas deste tipo. 

Número de pás — 60, 

As velocidades à entrada são completamente defi- 


nidas à custa destes valores, assim: 
Velocidades 


wu = 7,85 m/5 
= = 3,98 m/s 
w, — 8,80 m/s 


Algumas das velocidades podem ainda ser conhe- 
cidas à saída 


us = 19,62 m/s 


Wo = 1,59 m/s 


O ângulo da pá à entrada é 
B, = 26,98º 
cálculo são 


Outras constantes necessárias do 


ainda 


w=157 rad/s 
k = 0,199 m'/s 


Os resultados são apresentados sob a forma de grá- 
ficos nas figuras 13 a 24. 


10 — ESCOAMENTOS COM CONDIÇÕES DE ENTRADA 
E SAÍDA IMPOSTAS 


A — Pã circular 


Consideremos uma roda centrifuga com pás circula- 
res, figura 11. Designemos por O o centro da secção 
recta da roda, C o centro de curvatura da secção da pá, 
e P, um ponto sobre a pá à entrada. Os pontos O, Ce P, 
definem um triângulo em que um dos ângulos internos 


39 


5 é suplemento do ângulo formado pela velocidade re- 


lativa w, à entrada com a direcção tangencial. 
Os valores 
R— raio de curvatura da secção da pá 


| — distância entre os pontos O e € 

podem ser determinados se forem conhecidos os triân- 
gulos das velocidades à entrada e à saída da roda. 
Assim, pelo teorema de Carnot; 


EEE EPE 
= V Re+r—2 Rr cos B' 


(10.1-A) 


Consideremos o triângulo OCP, em que P; é um 
ponto da pá na periferia da roda, Conhecendo o triân- 
gulo das velocidades à saída o ângulo A'; fica definido, 
a então 


|=4/W+Hr;—2Rrncos f' 
Ta PD mam 


Igualando as expressões (10.1-A) e (10.1-B) ob- 
temos 


d 


Pi — Ta 
R>[-DD DDD > >—>———— 
2 (rs cos B':— ri cos 8º) (10.2) 
e a substituição deste valor em (10.1) permite calcular | 
Por aplicação sucessiva do teorema de Carnot ao 
triângulo OCP podemos também obter as relações (fi- 
qura 11): 


RE+r— po: 
RS e “2Rr ) (10.3) 
R+E—r 
nn (10.4) 
R+ Er 
PRSD (10.5) 
com 
ESabe 
CS ad to (10.6) 


À igualdade (10.5) é a equação das linhas de cor- 
rente. 
O perimetro das linhas de corrente calcula-se facil- 
mente pela expressão 
s = Re (10.7) 
Para o cálculo da velocidade e aceleração é neces- 


sário conhecer o valor do ângulo 5, e o valor da deri- 
vada df/dr obtida por diferenciação da equação (10.3) 


dA cosec à cotg 8 


— — Ex 
dr 2 R r (10.8) 


À pressão exercida sobre as pás calcula-se também, 
substituindo os valores agora obtidos na relação (4.16). 


B — Variação Linear de f 


Admitindo um crescimento linear de 5 ao longo da 
pá, podemos escrever, para o perímetro elementar da 
pá ds 


dB 
- = a = Constante 
dr (10.9) 


Então, por (2.30) 


dB 
———— = à cosec ;! 
dr (10.10) 
a integrando 
cosd)— cosg=-a lr-ri) (10.171) 


Esta relação permite-nos ainda calcular o valor de a 


cos B — cos 8; 
ã— ——— oia 
Fa — E (10.12) 
Não é possível obter a equação das linhas de co-- 
rente na sua forma integrada. No entanto, podemos fazer 
o seu cálculo aproximadamente 
AB cotg É 


nã. 


AT r (10.13) 
para intervalos pequenos À r. 


Substituindo (10.10) nas equações das velocidades 
a acelerações podemos ainda escrever: 


a 


Wo 


w' = — (nr+dr cosec B 81) 
r (10.14) 


c' = — lar cosecf w — uw—2uw )6 


(10.15) 
À diferença de pressão entra as pás será: 


2% 


“—— senf cos8 w (1 -—ar secs 
z 


- send + 2 t 
m a AM (1016) 
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Termodinâmica molecular do gás perfeito 


RESUMO 


Apresenta-se uma revisão critica da determinação 
teórica de propriedades termodinâmicas de gases per- 
feitos. Com base na Termodinâmica Estatistica e em 
dados fornecidos pela análise espectroscópica indicam- 
-se métodos de cálculo de várias contribuições para os 
valores das funções termodinâmicas. Finalmente, aplica- 
-se o estudo às moléculas de CO, NH, e SO, e compa- 
ram-se resultados calculados de S e € p com valores 
experimentais. 


1 — INTRODUÇÃO 


Às propriedades termodinâmicas de uma substância 
— sólida ou fluida —- estão intimamente relacionadas 
com as qualidades inerentes às moléculas que a consti- 
tuem. Por outras palavras, e , mais abrevidamente, pode 
afirmar-se que as propriedades microscópicas (das molé- 
culas, ou dos «sistemash) determinam os valores das 
grandezas macroscópicas (do conjunto, ou da «assem- 
bleiay). 


A asserção anterior, que é, aliás, intuitiva, levou 
a pensar na possibilidade do cálculo de propriedades 
termodinêmicas de uma assembleia, conhecidos os afri- 
butos que caracterizam os sistemas que a compõem. 
Porém, para passar do microscópico ao macroscópico 
é necessário dispôr de uma teoria que estabeleça a liga- 
ção entre os dois domínios. Ora a Termodinâmica Esta- 
tistica, na fase actual de compreensão das interacções 
moleculares, desempenha cabalmente tal função para 
o estado gasoso. 

Nos parágrafos que se seguem vamos ilustrar este 
ponto, através da determinação de calores específicos 
a volume constante, C, , e de entropias, S, de gases 
constituídos por moléculas simples, mediante a aplica- 
ção de vários modelos físicos e a utilização de dados 
moleculares, fornecidos pela análise espectrescópica. 


J. à. MARTINHO SIMÕES * e 
C. &. NIETO DE CASTRO 
Laboratório de Quimica-Fisica 


Instituto Superior Técnico 


ABSTRACT 


A critical review about the theoretical determina- 
tion of thermodynamic properties of perfect gases is 
made. Based in Statistical thermodynamics and in data 
derived from spectroscopic analysis, methods of cal- 
culus of several contributions for the values of the 
thermodynamic functions are indicated. Finally, the 
study is applied to the calculation of molar entropy and 
heat capacity of simple polyatomic molecules, viz CO,, 
NH, and SO. A comparison with exprerimental values 
is made. 


O princípio do método é simples. A termodinâmica 
Estatística introduz uma função denominada «função de 
partição», Z, a qual está relacionada, quer com os níveis 
energéticos s; das moléculas, quer com as proprieda- 
des termodinâmicas da assembleia (como o seu nome 
indica, Z determina a distribuição dos sistemas pelos 
níveis £ |. Assim, pode demonstrar-se que [1]: 


-W — go SK 
Z= Lwe e;/kT (1) 
i 
a que, por exemplo, 
E dInzZ 
= Tl mm (2) 
NkT dT VN 
A 
e E mb (3) 
NkT 
C, | o dln Z ; 
a di E Te ) 
Nk | 9T aT VN V, N 
(4) 


sendo E e À respectivamente a energia interna e a ener- 
gia de Helmholtz do gás, w, a degenerescência do nível 
*;, 1 a temperatura absoluta, N o número de Avogadro, k 
a constante de Boltzmann e V o volume da assembleia, 


“* Os resultados agora apresentados tiveram origem num trabalho realizado no ano lectivo de 1973/74, integrado na cadeira 


de Laboratórios, por J. A. M. 5. e M. Isabel 5. MARTINS. 
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Todo o problema parece, pois, resumir-se ao desen- 
volvimento de Z em função da temperatura, desde que 
sejam conhecidos os valores de «, e de w,. 

A, energia total, «., de uma molécula pode ser obtida 
pela resolução da equação de Schrôódinger 


Hy = ey (5), 


onde H representa o hamiltoniano. Uma maneira usual 
do simplificar este problema consiste em supôr inde- 
pendentes as várias formas de energia molecular: trans- 
lacional (e, ), rotacional (e, ), vibracional (=, ) e electró- 
nica (e, |. Assim, o hamiltoniano da equação (5) aparece 
como a soma de quatro termos independentes, o mesmo 
acontecendo, por consequência a e: 


EE ep e tt (6). 


A confirmação desta hipótese é feita «a posterioriy 
pela comparação de valores calculados com resultados 
experimentais. 

Os três primeiros termos de (6) correspondem, em 
geral, a outras tantas soluções da equação de 5Schrôdin- 
ger aplicada respectivamente ao modelo de uma partícula 
pontual de massa m, contida numa caixa de volume V 
(cue corresponde ao modelo do gás perfeito), ao modelo 
do rótor rigido e ao modelo do oscilador harmónico 
linear. À resolução de cada um destes problemas conduz 
a uma expressão quântica dos níveis de energia. 

A equação (6) pode ser sofisticada. Por um lado, 
6 possível introduzir um termo que dá conta da inte- 
racção entre a rotação e a vibração da molécula (ey! 
c que equivale a admitir a interdependência entre estas 
duas formzs de energia. Simultaneamente, qualquer um 
dos três modelos físicos referidos pode ser refinado, 
medianta a introdução de correcções, como as devidas 
à anarmonicidade (do oscilador), e,, à distorção centri- 
fuga (do rótor), egçr StC.. Vem então 


= -+- E + Es + Ea + EG Ea Ega rios (7). 


Introduzindo (7) em (71) tira-se que a função de 
da partição será o produto de sete (ou mais) termos, 
cada um correspondente a um «tipo» de energia. O mes- 
mo se poderá dizer a respeito das equações (2), (3) e 
(4), onde os valores de E, A e €C, virão dados por uma 
soma de sete parcelas. 

Em conclusão, conhecidas as exoressões quantifi- 
cadas dos termos de (7) e desenvolvendo os somatórios 
de (1) é possível o cálculo das contribuições das várias 
formas de enregia molecular para as propriedades termo- 
dinâmicas do gás. 


| 


2 -— ESTUDO DAS PROPRIEDADES 
TERMODINÂMICAS DE MOLÉCULAS 
SIMPLES 


Z:1 


Como foi referido na introdução, a Termodinâmica 
Estatistica indica-nos a forma como os sistemas se dis- 


1-— O (x) significa «função da ordem de x». 
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tribuem pelos diferentes níveis de energia acessíveis 
O número de moléculas n; num nivel E; será 


w 6 — E fkT 
n;=N. 5 (8) 

Quando se faz incidir sobre uma amostra de gás 
um feixe de energia «, as moléculas são excitadas a 
níveis energéticos superiores. Como se sabe, a espectros- 
copia do Raman e infravermelhos permite-nos dar conta 
de dois tipos de niveis de energia — rotacionais e vibra- 
cicnais — e a do visivel e ultravioleta mostra-nos a exis- 
tência dos niveis electrônicos. 

A energia necessária para «promover» as molé- 
culas do estado fundamental ao primeiro estado exci- 
tado, bem como as diferenças entre os niveis corres- 
pondentes às várias formas de energia, varia do se- 
guinte modo: 


Ar; SS dr, «SS de, << dr, 


Na ausência do radiações incidentes (gás em equi- 
librio) as moléculas só poderão ser promovidas pela 
egitação térmica, cuja energia é da ordem de kT. Sendo 
assim, a análise qualitativa das diferentes contribuições 
para es propriedades termodinâmicas poderá ser feita 
através da comparação do valor de kT e do valor de As 
correspondente aos graus de liberdade em causa, Então, 
por exemplo, nara uma molécula diatómica típica, tem-se, 


Ac, =0 (10 ev) 

Ae, =0 (5x 10 *ev] 

Ae = O (0,3 eM) 

de, =0 (5 eV) 

kT (300 K) = 0 (0,03 eV) 


ua 


Fig. 1 


sabendo que dr = ejry — e; (Figura 1), pode- 
mos tirar algumas conclusões com base na relação 
(8). Assim, por exemplo à temperatura ambiente, cons- 
tatamos que os níveis translacionais estão todos ocupa- 
dos (As, << kT, de onde a aproximação de considerar 
um econtinuo» de níveis de energia e a igualdade dos 
resultados obtidos pelas teorias clássica e quântica); os 
níveis rotacionais excitados estarão também bastante 
ocupados (Ae, é cerca de 60 vezes menor que kT); 
os vibracionais superiores estarão praticamente vazios 
(As, é cerca de 10 vezes superior a kT); finalmente, 
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no caso dos niveis electrónicos verifica-se apenas a 
ocupação do estado fundamental, Es: Portanto, se para 
uma dada temperatura compararmos A: com kT pode- 
remos ter ideia das diferentes contribuições para as pro- 
prisgdades termodinâmicas, ou seja, mais precisamente, 
se será válida a aplicação do Princípio da Equipartição 
da Energia, o qual, como se sabe, é apenas utilizável 
no caso clássico (Ae << kT) [4]. 


2.2 — CONTRIBUIÇÃO TRANSLACIONAL 
A resolução da equação de Schrôdinger aplicada 


ao modelo de uma partícula pontual de massa m que 
se move num volume V conduz a [2] 


h? p* q Pê 1 
“= = +t—— +— | (9) 


(10). 


tal que abc = WV.h é a constante de Planck e p, q, r são 
números quânticos que tomam valores inteiros, positivos, 
superiores ou iguais à unidade. 

A substituição de (9) em (1) conduz à relação 


Do o DO 
| a pº q 
| A 


D=4 0=1 1=79 à 


rá hº 
+ cº hair! 


a qual pode ssr factorizada num produto de três soma- 
tórios, em p, em g e em r. À resolução de cada um dos 
termos é imediata se nos lembrarmos de que de, << kT, 
mesmo a temperaturas apreciavelmente baixas. Nestas 
condições, é fácil ver que as quantidades 


(11), 


hº 
Ca = Bm a” kT 
hº 
“b= gm b? kT 
hº 


“e= 8&m c kT 


são muito inferiores à unidade e que, portanto, será lícito 
substituirmos os somatórios por integrais [1]. Vem 


oo ai =] =| 
Z | == [ eaptdo | E — ap da. | e a dr 
“a a PT 
(12), 
o que nos conduz a 


(13), 


quo é, pois, a expressão da função de partição transla- 
cional para uma molécula de um gás perfeito, 
A obtenção das expressões das propriedades termo- 
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dinâmicas é imediata. Contudo, não podemos esquecer 
quo estamos a estudar uma assembleia composta por 
sistemas idênticos, independentes e não-localizados. Ora 
demonstra-se [3] que, neste caso, a função de partição 
translacional da assembleia dos N sistemas é 


ZN 
o,— ss (14). 
N! 


À única alteração a fazer é na equação (3), a qual 
se pode escrever sob a forma geral 


A=-kT In Z2N 


=-—-kT In Q (15). 
Assim, para O caso em estudo, 
A=—kT In ZN+kT In NI (16). 


Como N é muito grande (— 10%) é válida a apro- 
mação de Stirling [1] 


In NI=N In N—N (17). 


Substituindo-se em (16) e rearranjando, obtém-se 


A 
——) =-nZ, +InN-1 
NkT |), 


em que Z, é dada por (13). 


(18), 


A relação final de (A/NkT), aparece numa forma 
mais prática de utilização se recorrermos à equação de 
estado do gás perfeito, pV = NkT, e, simultaneamente, 
se efectuarmos algumas transformações, de modo a 
podermos substituir as várias grandezas, que nela figu- 
ram, em unidades mais correntes. Vem então 


A 3 3 5 T 
=j=""=>"nM=-— n—-h'p 
NkT | 2 2 2 1,594 

(19), 


sendo M o peso molecular (em gramas), T a temperatura 
absoluta e p a pressão (em atmosferas). 

às equações (2) e (4), referentes respectivamente 
a EeacC, , não necessitam de alterações, uma vez que 
envolvem derivadas de In Z, em ordem a T (o que anula 
o termo In Nl!). Por consequência, substituindo (13) em 
(2) o em (4) e resolvendo as derivadas tira-se que 


(É 3 
Mem 
“NkT 


(21 IP 


resultados que são concordantes com os obtidos com 
base no Princípio da Equipartição da Energia (três graus 


1 
de liberdade vezes -——KkT para a energia de cada um). 
2 


2.3 — CONTRIBUIÇÃO ELECTRÔNICA 


Ão contrário do que acontece para o caso da trans- 
lação, agora de. >> kT, mesmo a temperaturas eleva- 
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das. Vejamos quais as consequências deste facto. Escre- 
vendo alguns termos de (1), vem: 


É. = eokT + w é E [kT E grid (22). 


É 


Como, por convenção, c,, = O, segue-se que Ag = 
= e >> kT. Assim, o segundo termo da equação ante- 
rior é praticamente igual a zero. E, por maioria de razão, 
O mesmo acontece aos seguintes (em =, e, etc.), cuja 
diferença em relação a «, é ainda maior. O resultado 
final é, pois, 


e 


ZW (23): 


Com o auxílio de (2), (3) e (4) obtêm-se as contri- 
buições electrónicas para as propriedades termodinâmi- 
cas respectivas: 


E 
) = O (24), 
NkT /. 
A 
= À 25 
NkT ). fic ER 
C, 
= (26). 
Nk É 


2.4 — CONTRIBUIÇÃO VIBRACIONAL 


A contribuição vibracional é deduzida, em primeira 
aproximação, com base no modelo do oscilador harmó- 
nico linear. 

A resolução da equação da Schrôringer aplicada a 
esto problema conduz à expressão: 

“/ 


à dO 
Ev: = (mu E E 2 ) ho, (27), 
1 v . 


onda o somatório se estenda ao número de frequências 
1; da molécula”; d; é a degenerescência da vibração 
RE: N, é o número quântico da mesma vibração, o qual 
pode tomar valores inteiros, positivos, superiores ou 
iguais a zero. 

A substituição de (27) em (1), após tratamento 
adequado, leva-nos à expressão de Z, [1]: 


= ai (28), 
' 11-67 hr;/kT ) 


dao onde se tiram as contribuições para as propriedades 
termodinâmicas; 


E o “ hr, | hr;/kT E (29) 
NkT /, Tlar "e hykT fi ' 


hr. 
em (1-0 Me /6T) dy 
i 


2kT 


(30). 


2.5 — CONTRIBUIÇÃO ROTACIONAL 


O estudo da contribuição rotacional tem que ser rea- 
li ado separademente para vários tipos de moléculas, 
cada um deles caracterizado vor um determinado com- 
portamento de rotação. Nomeadamente, a subdivisão é 
feita atendendo aos valores dos momentos de inércia 
(1, ; lho le) a sua inter-relação. 

Vamos apenas abordar o problema de moléculas 
lineares (di- ou poliatómicas), de moléculas de cuma si- 
métrico o de moléculas de cume assimétrico. Para estes 
três casos tem-se, respectivamente: 

Hace =0; |, + O 


é] ) 


= == H | E 
A Er lei Lo po, Cc 


2.5.1 — MOLÉCULAS LINEARES 


Bb, expressão dos niveis de energia do rótor rígido, 
obtida pela resolução da equação de Schródinger, é, 
para este caso, 


= Bh Jd (+) (32), 


sendo a degenerescência do nível E] 
A = 2)+ 1 (33). 


J é o número quântico rotacional, que toma valores intei- 
ros, positivos, superiores ou iguais a zero; B, é a cha- 
mada constante rotacional para o estado fundamental 
de vibração (y = o), dada por 


h 


B, e ag ar? | (34), 


tu 


onde |, é o momento de inércia correspondente a y— o. 
A função do partição tem, pois, a forma inicial: 


DO 


Z Ss (us + 1; e Bah + WIKI (35). 


r cer 
J- O 


Chamemos u à razão 
B,h 
kT 


a escrevamos es primeiros termos do somatório (35): 


u= (36) 


Z=1+3072+5076U+... (37). 


Se u>1 (ou T<B, h/k) a sério converge com 
rapidez, a qual será tanto mais acentuada quanto menor 


(2) Para cuaçÃos om mepouse, infinitamente afastados, Z , =w, exp (D /kT), sendo D, a energia da dissociação espectroscópica, 
(3) dama ro res utilizam o somatório estendido aos modos normais da vibração da molécula (3n-5 para moláculas lineares e 3n-6 
para as fr Mtantes, sendo n o número de átomos), pelo que, neste caso, será d = 1 para todo o |. 
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for a temperatura, Portanto, é possivel desprezar termos 
superiores a uma dada ordem e, em seguida, deduzir as 
expressões das propriedades termodinâmicas. No en- 
tanto, a razão B, h/k é, em geral, bastante inferior a 
10 K. O estudo deste caso não tem, assim, qualquer inte- 
resse imediato, uma vez que, a essas temperaturas, a hi- 
pótese do gás perfeito perde a sua validade. 

Seu< 1 (ouT>B, h/k) a série converge muito 
mais lentamente que no caso anterior. Por isso não é 
viável trabalharmos com a expressão de AR tal como se 
apresenta em (35), dado que, para tal, seria necessário 
considerarmos um número enorme de termos (tanto 
maior quanto maior a temperatura e o grau de precisão 
desejado), o que complicaria extraordinariamente as re- 
lações das propriedades termodinâmicas, 


à maneira de resolver o problema consiste em apli- 
car à equação (35) o teorema de Euler-MacLaurin: 


Do 1 DO 
o» f (x) — | 


x=0 : 


ft (x) dx+n (38), 


onde q representa o valor do erro que se comete quando 
se substitui o somatório pelo integral, e que é dado por 


1 1 
ii = f = 2) E! g (iii) (0) — 
] 2 (o) 13 to) | (0) 
1 
e f Iv) toj+... (39). 
3040 
Fazendo, no caso em estudo, 
fix)=(2x+1) em xtx+1)u (40), 


e recorrendo ao teorema citado, obtém-se facilmente que 


1 u Ti 4 u” 
+ (1+ 3 "45 * 31 
u! 
E oc) (41), 


desprezando termos de ordem superior. 


Notemos que, quanto maior for a temperatura, me- 
nor será o valor de u e, portanto, maior o número de 
termos de (41) que podemos desprezar (recordemos 
que estamos a analizar o caso em que u é sempre in- 
ferior à unidade). No caso limite, isto é com u << 1, 


(42), 


que é a função de partição rotacional clássica. 

Em resumo, os termos polinominais de (41) devem 
ser utilizados quando a temperatura não é suficiente- 
mente alta para que, com um dado grau de precisão, se 
possam obter valores de Es 


(4) 


O passo que vamos dar em seguida é, uma vez 
mais, o cálculo das expressões termodinâmicas. 

No que se refere à energia de Helmholtz não existe 
qualquer problema. De acordo com (3) tem-se 


A u u” 
fm ) = Ino + Inu In (1+ L + 
NkT F, 3 15 

4 u” u 
— fe | 43). 
E 315 315 ) 8) 


Para a energia interna e para o calor específico, 
ondo se torna necessário o cálculo da derivadas, pode 
pôr-se uma questão quanto ao procedimento a adoptar. 
Por um lado, podemos desenvolver o In Z, em série de 
MacLaurin, desprezando os termos de ordem superior 
à quarta” e, em seguida, derivar; por outro lado, é pos- 
sível derivar directamente o In Z, . Embora a via que 
pareca mais correcta seja esta última, o que se verifica 
é que, no primeiro caso, ocorre uma compensação de 
erros (devidos às aproximações de (38), (41) e ao de- 
senvolvimento de MacLaurin), obtendo-se, vor conse- 
quência, um valor mais correcto. No entanto, e como 
seria da esperar, registemos que a diferença entre as 
duas opções só se manifesta a temperaturas razoavel- 
mente baixas (— BO K). 

Nestas condições, efectuando-se o desenvolvimento 
somos conduzidos a 


u u 8u' 59 uí 
In Z=-—-Ine-— In vu+ — + 4———— + 
i 3 90 2835 56700 
(44). 
Assim, recorrendo a [2) e a (4), vem 
E u u” 8 u 59 u' 
| | Su esse ES (45) 
NkT !, 3 45 945 14175 
a 
f My ué 16 u 59 u' 
| ) =1+—0 A (46). 
Nk 45 945 4725 


E E 
No caso deu <<1,(C, /Nk), =1e (E/NkT), 
= 1, o que coincide com o resultado clássico, seja o pre- 
visto pela aplicação da Mecânica Clássica (equação 
(42)), seja obtido através do Princípio da Equipartição 
da Energia. 


2.5.2 — MOLÉCULAS DE CUME SIMÉTRICO 


A relação dos niveis energéticos do rotor rigido é 
agora 


e =B,hJN + +HhIA,—B)K (47), 


sendo A, a constante rotacional correspondente a a 
para r= O o K um número quântico que toma valores 
inteiros compreendidos entre J e-J. SeIKI<SO os 
estados são, portanto, duplamente degenerados. 


é o chamado «número de simetriar [1], que é, neste caso, introduzido «ad hoc» para dar conta da simetria de rotação da 


molécula: o= 2 para moléculas diatômicas homonuclgares ou moléculas poliatómicas com um plano de simetria norma! ao eixo inter- 
nuclear e que passa pelo centro de massa; q = 1 para os casos restantes. Para os tipos de moléculas de que trataremos nos parágrafos 


seguintes, q pode tomar valores superiores a 2. 


(5) De modo que haja coerência com o facto de, na equação de Ae termos desprezado os termos de ordem superior à quarta. 
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ql 


O tratamento realizado neste caso é semelhante ao 
adoptado para as moléculas lineares. Assim, a expressão 
de £ que se obtém, e que foi pela primeira vez dedu- 
zido por Viney [5,6] é 


1 Tom skT ST 1 
2, eme WE gd dt m— 
a B. Ao h E 12 
7 N 
au + -— "e df cg) | (48), 
480 
onde u é dado por (36) e 
B 
a=1-— (49). 
A, 


As opções que se podem pôr, com vista à determi- 
nação das propriedades termodinâmicas são as mesmas 
que as referidas no 8 2.5.1. Contudo, o que agora se ve- 
rifica é que a derivada directa do In Z, conduz a re- 
sultados ligeiramente melhores, embora a diferença entre 
os valores obtidos pelas duas vias seja, neste caso, 
ainda mais desprezável, mesmo a temperaturas bastante 
baixas. Assim, as expressões a que se chega são 


E u 3 
NKkT É a 
1 7 
-— e 0 e] a u” 
12 480 
q E 50), 
1 7 (50) 
7 + ——— ay ty + a u 
12 480 
A u 1 
- ) In o PR ip 
NkT te 4 2 
VT kT 45 1 7 
In e) (——) Un t— uy 5 o u*) 
BA, h 12 
8 (51) 
Ga 3 
[aÃ sas 
Nk , 2 
1 7 49 
—— q = e] atas in a“ ut 
45 2880 115200 
logs 
1 13 7 49 
1+t—qut-—— eu+——— qgu+—u' 
6 360 2880 230400 
(52). 


Notemos finalemnte que, se u << 1, somos condu- 
zidos aos resultados clássicos de (E/NKkT) rede (C, 
[Nk), . 


2.5.3 — MOLÉCULAS DE CUME ASSIMÉTRICO 


Ao contrário do que se passa nos tipos de molé- 
culas que estudámos até aqui, não existe, neste caso, 
solução da equação de Schrôdinger aplicada ao modelo 
do rótor rígido. 

A dificuldade pode, no entanto, ser ultrapassada, 
mediante a adopção de uma aproximação, ou, ainda, pela 
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aplicação de métodos dedutivos da função de partição 
que não requeiram o conhecimento explícito dos valores 
próprios do hamiltoniano. 

A aproximação mencionada consiste em tratar o 
problema de uma molécula de cume assimétrico como 
se esta se tratasse de uma molécula de cume simétrico. 
Nomeadamente, atendendo a que, em muitas moléculas 
daquele tipo, B, — C, . as expressões de Z, e das 
funções termodinâmicas serão idênticas às obtidas no 
parágrafo anterior, tendo apenas como única diferença, 
o facto de os valores de u e de « virem agora 


= N B, Co h (53) 
U— 
kT 
| BE 
o q NV Toto (54); 
xx == = 
f 


Q 


isto é, substitui-se B, pela média geométrica entre 
Do O GC - 


Todavia, existem moléculas para as quais a aproxi- 
mação descrita se revela bastante grosseira, na medida 
em que B, é substancialmente diferente de C, . Se 
assim for é necessário enveredar pela segunda opção. 

Stripp e Kirkwood [7] deduziram uma expressão de 


Z, avropriada para moléculas de cume assimétrico, a 
qual! se reduz à equação (48) quando B, —C, 
1 q kT (5 712 
Z. —— ) 
o A Bolo h 
1 1 
ap as RE = yt) 55). 
72 *** 280 | Epsa 


onde C, é a constante rotacional correspondente a |, 
para v=0O, 


A,B, B,C, CA 
por . E] ú o t q Ú 
p=2A,+20,+ 26 a 
o) E) o 
(56), 


y = 10 (A) FB +C, ): 12(A,85+ B,Cot CA, 


| A o B, TABS | Bo C,+ Bo Co 
C,, A, 
| Co Ao Tt CAS 43 I A q Bo Ma Qi | coro | 
B,, | | C, Pia B, 
(57) 
BG 
Rr (58). 


Com base em (55) podemos, pois, obter as expres- 
soes termodinâmicas. Tal como no caso precedente, os 
cálculos por nós efectuados também aqui revelaram li- 
geiramente favorável a via da derivada directa do In 


AR - Vem então 
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1 1 E 
-—- Bpt-— ye =2r7+ ET PT” ' 
E 3 12 B 240 Y Ee NKT In r 20r T'+ 296 r'T (65) 
RD eai a OI, 
fra, 2 1 1 au: A: 296 
+—— ptt— y — =-2+T-10"T'-— PT' 
2 "tao" ER E E RR 166) 
a 
A 1 q KT: 4 dr E 
— | —= In PE ——— | l ) Ss PR | 3 a rã 
fd] lr Fera E | NO) 4 4T+ 60 PT + 1184 PT (67). 
1 1 (60) 
= tm (1 Fa + zo 1º ) 
à 2.6.2 — MOLÉCULAS DE CUME SIMÉTRICO 
B 
Cc 3 Para moléculas de cume simétrico 
o) 
a) pi s,=BohJ (J+1)+(A,— B)hK — 
1 1 1 =D bh (+18) -D hJ J+HI)K?-DEhKkK' 
| = Ye J | JK | | k 
ago 1º 2880 Br + 715200 115200 (68), 
1 + E fe as )e+ pyt E representando D | é De Dy as constan- 
240! “Tas 2880 230400 (61). tes de distorção centrífuga. 


Como é fácil verificar, para temperaturas suficien- 
temente altas obtêm-se os valores clássicos de (E/NkT) ; 
e de (C, /Nk), 


2.6 — CONTRIBUIÇÃO DA DISTORÇÃO CENTRIFUGA 


O método de obtenção das expressões termodinã- 
micas para este caso é análogo ao que estudâmos no 
5 2.5. Apenas há a referir ligeiras alterações de carácter 
matemático, que não valerá a pena discutir, até porque 
facilmente se podem encontrar na literatura indicada. 


2.6.1 — MOLÉCULAS LINEARES 


Os níveis de energia do rótor não-rígido são dados 
por 
e-= BhJ + 1]= Db” (J + 1) (62), 


sendo D, a constante de distorção centrífuga para 
v = O. A degenerescência de cada nivel rotacional con- 
tinua a ser igual a 2J + 1. 

McDowell [8], baseado no último termo do se- 
gundo membro de (62) (que dá conta da distorção cen- 
trifuga) e aplicando o tratamento referido, chegou a 


In Zgç= 1 + 2rT + 12r Tº+ 120r T' (63), 
onda 
na k Do 
hB, (64). 


Os termos omissos na equação (63) são despre- 
záveis mesmo a temperaturas altas pois o valor de r é 


bastante pequeno (— 10 KT!) 


As propriedades termodinâmicas foram também 
calculadas segundo o processo usual (neste caso, de- 
senvolvendo em séria e derivando), chegando-se a 
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Igualmente para este caso McDowell [8] chegou 
às equações de In Z,. e das propriedades termodinã- 
micas, tendo como ponto de partida os três últimos 
termos do segundo membro da equação (68): 


1 | 
nZaç= e, T+ 2p=e3)r (69), 
f E 4 
es =, T+ — ET 70), 
UN E É bs) di 
A q Rr 
—— |) =-4T—-— [ 2m-—pº) Tº 7 
ae. É E da e) T dd: 
a 
C, j 
=207M+3 (2 »—- te 2), 

Nk IN pa pa 04] (72) 
sendo 

: “(8+4 +3m') D;+ |2 
Pp =—] m m - (2m+4 
"ans; : 

+ 3m* ) D+ 20º Dx | (73), 

3 k? 
E = —— sor | (128 + 64m + 48 mº* + 40m' + 

32h'B, 


+35 m') 03+2 (16m + 24m* + 30m' + 35 m?) 


1 
+ Djxl +2 (10m+ 35m') D gDy + 35 m'Dg | 
(74) 
B 
o m=— (75). 
As 
73 


2.6.3 — MOLÉCULAS DE CUME ASSIMÉTRICO 


Como vimos no 4 2.5.3, a contribuição rotacional 
para as funções termodinâmicas de moléculas de cume 
assimétrico não pode ser estudada pelos métodos usuais 
(uma vez que não existe solução da equação de Schrô- 
dinger aplicada ao modelo do rótor rígido), passando-se 
o mesmo com respeito ao problema da distorção cen- 
trifuga. 

Uma maneira de resolver esta questão, ou seja, de 
calcular o efeito da distorção centrifuga nas proprieda- 
des termodinâmicas das moléculas em estudo, será a de 
adoptar um tratamento baseado na via clássica. Con- 
cretamente, o assunto pode ser abordado encarando o 
problema do rótor não-rígido através da Mecânica Clás- 
sica. 

Esta aproximação foi realizada por Wilson [10], 
que assim conseguiu obter uma expressão de E dc 


Zac 1 +U'T (76) 


onde É' é uma função das constantes de força da molé- 
cula ( 7) a tem normalmente um valor da ordem de 
10-5 K—1: 


1 | 
E À E 
= k( 3 e AG + 3 ogyy Bo TSC razao! 


+ 2Cexyy A, B, ai 2Coyzz Bo C, Ei 2 E rena Cc, Aç+ 
Eds Cxvyxy A, Bo + 4z vayzBo Cs + SE ue Co Ao) 


(77). 


Da (76) Wilson tirou as expressões termodinâmicas 
desenvolvendo o In Zac em série a derivando. Chegou a 


E 
(57) =" (18), 
NkT Ju 

A 

DA, age 9 
NkT ac 

, C, 

6 a )=2e T (80). 


2.7 — CONTRIBUIÇÃO DA ANARMONICIDADE 


Os níveis de energia do oscilador anarmónico li- 
near são dados por 


d. a d. 
. pe é Ene E 
ed (y "a hr; na(s, "3 X hr; + 
d. d 
2, à ane E [O ias 2 
nm (v a )( É + 2 x;h v; Fe. 
(81), 
ondo os dois últimos termos do segundo membro, atra- 
vês das chamadas constantes de anarmonicidade, Xj 
OX. dão conta deste fenómeno, Tal como no & 2.4, 


os somatórios estendem-se ao número de frequências. 
O processo de tratar o problema, de modo a obter- 
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mos expressões de In E e das propriedades termodinã- 
micas é análogo ao utilizado nos parágrafos anteriores. 
*, diferença é que, neste czso, ele se revela muito mais 
trabalhoso. 

Woolley [11], com base em (81), chegou a uma 
expressão de In Z, onde apresenta, não só vários termos 
respeitantes às constantes x, e Xj - como também 
termos devidos às constantes de anarmonicidade de 
maior ordem (omissos na equação (81)). Na prática, 
alguns dos termos de Woolley têm uma contribuição 
desprezável, nelo que podemos eliminá-los. Por maioria 
de razão, as mencionadas constantes de ordem superior 
terão uma contribuição ainda menor, para além de que, 
p=ra a quase totalidade das moléculas, os seus valores 
não se encontram determinados. 

Pennington e Kobe [9] abordaram também ese 
problema, embora a expressão final a que chegaram 
contenha apenas os dois primeiros termos dados por 
Woolley (que são, no entanto, os de maior peso). 

A equação de In Z, que vamos utilizar é, pos, 
a do Woolley, reduzida aos termos de maior importância. 
Alterámos, porém, a sua forma, mediante o emprego de 
funções ly, definidas »or Pennington e Kobe, as quais 
nos facilitaram grandemente a obtenção das relações ter- 
modinâmicas e os próprios cálculos. Tem-se então: 


38 
E À a 


E 


” 
P; Pi + cm à detd; + 1) Pj 
i<j<k 
= | | Ez ] : 


l 1 1 
| A Fr AG SA | 
+ Xi Xi did; de “o; 0; Pk | + 


es [2%4%; d(d; +71) dy "q; "q; a) (82), 
Ra 


ER LA 2 
f; F; ( F + Fj ] | jolie [id (d, - 1) Ug; 


=te0| tm ta no +Xid, (d; + 1)ºg; 


1 1 
i E = Ê 
po +Hid+1) o] E 
à 1 
gp to, Spa Md E 
ROO TORA it 
1 1 
5 Ea ado Ma A a PRESA 
Ae vt Sd AL" 01 + E 


de nto dy (oe) | + 
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as sa add: | E e e TE Sem, Re Pe 


a 
7 “ 7 q dum 
Pk FP; Dj Pr Pr JRR E 

Eis 


=] 


| 2 X;Xud;(d +41) 


d, 'p; “gp; “q; 'p; | 'pilto; “pi ) | 
+ RP O) 9 | dia 
A 
| =-—InZ, (84) 
NkT |, a 
e] 
(1) 2 cit: RA SAO + Xdid 
E — oa pao md z E: ' é 
Nk fa e] at Muda Ra Mae 
“q; q; | ( “q; + “p; Ee 0 q; + “g) 1 | E 
“| ã no BR 
+ sá [Mi Srta AO pt (d+ 1); 
le <K : 


+ A 0) mo + (d; + 1)ºgq; 


— 


- Ra E E | : A 
F| Ff Pi pj | 2 d, E] pj P) P; + 


| | Ma fa 5 fa ay T; x 
+ Xj Xw8,4, 8, | Fi P) Per PP; pl pj + 


tal que 
X.: 
fa = — H (86) 
É 
X.s 
6 Ki = pas (87). 
H (n; pj)! po 


As funções Ny são definidas por 
g= 48" 17 
'p=s 0!(6"=—7)T* 
p=ste'(o8'+1) (868— 1)7? 


gp =2s (68-1)-* 


'p=28s (25e8-e*+1) (68—-1)* 
p=4s est(2s 0!—-2 ef+-s+2) (68-11) 
sC(es— 1) 


q=— (28 -et+1) (ui 


|| 


“m 


“mp — sº e* (e* — 1)? 


dg =3 (68— 1)! 


Ip=-In(1-06 78) 
ip=ste'(-20º +s 0 +4s 0%+2 
+s) (68— 4) (88), 
onde s é dado por 
hp 
“e CEE (89). 


2.8 — CONTRIBUIÇÃO DA INTERACÇÃO 
ROTAÇÃO-VIBRAÇÃO 


Sao uma molécula passa do estado fundamental de 
vibração para um estado excitado, a distância média 
entra os seus átomos varia. Mais precisamente, se Ea 
representar a referida distância paravV =O e Fr, as 
reportar a um nível energético mais elevado, tem-se que 
t, >, - Sendo assim, para os momentos de inércia 
verificar-se-ã a relação |, > |; ; e, para as constantes 
rotacionais, será, por exemplo, B. e 

É precisamente neste fenómeno que consiste a in- 
teracção rotação-vibração. À variação de B com o nú- 
mero quântico vibracional pode traduzir-se pela expres- 
são 


6? 
B,= B, (1 SE Do, 4 ] (90), 


estendendo-se o somatório ao número de frequências 
da molécula e sendo b, dado por 


= + (91); 


et E são as chamadas constantes de interacção rotação- 


vibração”. 

Introduzindo as equações do tipo de (90) — para 
A,. B, o C, — nas expressões do rótor rígido e atra- 
vés de um tratamento adequado, Pennington e Kobe [9] 
lograram obter expressões da contribuição em estudo 
para vários tipos de moléculas. Chegaram a 


(92), 


onde o número de termos do somatório iguala o número 
de frequências F; é função de a, «b, ec, ea relação 
que o traduz varia consoante o comportamento rotacio- 
nal da molécula. Assim, para moléculas lineares, 


16) Para as constantes rotacionais À e C existem relações idênticas a (90 e a (91), definindo-se as grandezas aj = ab 7a, 


e C; a af 1Co 
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To 


fr= b + by (93): 
para moléculas de cume simétrico 
a; b; 3a? 


I j 
+ f (94); 
2 8 


Mis: PE b? 
F; b, + 2 + D; + 


para moléculas de cume assimétrico 

a+bi+te a +bitcç at b,+ cj) 

f; = — e - em am 
2 4 8 

(95). 


As contribuições para as propriedades termodinã- 
micas vêm então dadas por 


E "| 
| ddr, toi (96), 
NkT Ju 
A “ 
aee] =— dr; (97) 
 NkT de 
| 6, “ 
o [e ) = dir q; (98). 


d; é, como vimos, a degenerescôncia da frequência 
vw; as funções " q são definidas do modo apontado no 
parágrafo anterior. 


3, — CÁLCULOS 


Vamos agora ilustrar a teoria estudada através da 
determinação de calores específicos e de entropias mo- 
lares para as moléculas de CO, (linear), NH; (cume si- 
métrico) e SO, (cume assimétrico). 

Os cálculos foram realizados no computador do 
Centro de Cálculo da Universidade Técnica de Lisboa, 
após a elaboração de programas em Fortran IV para os 
casos apontados. Cada um destes programas permite 
obter separadamente todas as contribuições analisadas 
para as energias interna e de Helmholtz, para a entropia 
e para o calor específico a volume constante. 

Os valores selecionados para os cálculos são apre- 
sentados nas tabelas | e Il. Não houve possibilidade de 
obter dados das constantes de anarmonicidade do NH, e 
das constantes de interacção rotação-vibração do SO.. 
Por outro lado, a contribuição da distorção centrífuga 
para esta última molécula deve considerar-se como uma 
aproximação grosseira, uma vez que também desconhe- 
cemos o valor correcto de ?' (embora saibamos, no en- 


tanto, que ele é da ordem de 10 *K”!). 

Na tabela Ill mostram-se os valores das várias con- 
tribuições para o C, das três moléculas mencionadas, 
a algumas temperaturas. Na tabela IV comparam-se os 
valores calculados (somatórios de todas as contribui- 
ções) com resultados experimentais de C, ede Ss 
(a p= 1 atm). 


4. — CONCLUSÕES 


Como se pode observar na tabela |V, a concordân- 
cia entre os valores é muito satisfatória. Apenas para o 
NH,, à temperatura de 300 K, se manifesta um desvio 
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de cerca de 3 % em relação ao valor experimental; no 
entanto, para T = 500 K o erro baixa para 1 %. Assim, 
é lícito supor que a origem da discrepância reside no 
facto de termos considerado válido o modelo do gás per- 
feito, o que será, neste caso, pouco correcto, devido 
à possível presença de interacções moleculares signi- 
ficativas. 

Da tabela Ill pode inferir-se que a importância das 
contribuições devidas à anarmonicidade, à interacção ro- 
tação-vibração e à distorção centrífuga só se tornam sig- 
nificativas a altas temperaturas. Com o objectivo de 
darmos uma imagem mais clara do fenómeno, elaborá- 
mos os gráficos da figura 2, onde se mostra a variação 
com a temperatura de cada uma das contribuições para 
o valor global de C,. 

Pela análise das diferentes contribuições termodi- 
nâmicas verifica-se que todas elas são função da tem- 
peratura, à excepção de Z, , que é também função do 
volume, devido a estar intimamente relacionado com os 
graus de liberdade externos. As contribuições devidas 
a graus de liberdade internos são, assim, apenas função 
da temperatura e, evidentemente, das propriedades de 
cada molécula, 

Desejamos finalmente acentuar que o cálculo teó- 
rico de propriedades termodinâmicas com base em 
dados espectroscópicos pode ser realizado de um modo 
ainda mais sofisticado, ou seja, é possível entrar-se com 
um maior número de contribuições. É o caso, por exem- 
plo, das ressonâncias de Fermi e de Darling-Dennison 
[13], do efeito de anromonicidade na rotação, etc. 
Contudo, os seus valores são de tal modo desprezáveis 
que, para uma previsão de uma propriedade experimental, 
não merece a pena tê-los em atenção — para além de o 
cálculo ser relativamente complexo. 

Pretendeu-se dar uma visão sobre o cálculo de pro- 
priedades termodinâmicas do gás perfeito, entrando em 
linha de conta com os fenómenos moleculares mais re- 
levantes. Dada a actual acessibilidade do cálculo com- 
putacional e as modernas exigências do processo tecno- 
lógico e sua optimização, pareceu-nos da maior impor- 
tância desenvolver estes cálculos, fornecendo ao mesmo 
tempo um contributo para o estudo da termodinâmica 
molecular no |.S.T. 
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(al Valores de [12] 


TABELA | 


(b) Para o CO, valores de [13]; para o SO, valores de [14] 


(c) Valores de [5] 
(d) Valores de [15] 


(e) Para o CO? valores de [15]; para o NH. valores de [8] 


0.00121 
0.00072 
0.00309 . 


0,177 
0.015 
0,176 
0.230 


3 


tp — 
Tg 
LM] 


EE 


Vi /m 


44.011 


17.032 


64.06 


TABELA II 


(CN) HIC,),FIC ITANkK] (OC, /NKD,, | CC, /Nkdac | (C, /Nk)A 


2.5125 |— 23X 107 0.0002 


3.3262 — 0.0014 0.0007 
— 0.0022 
0.0003 


24X 10% 


TABELA IV 
S/cal mole”! K”! C,/cal mole K 

T/K | 
exp. (a) calc. exp. (a) calc. 
293.16 50.864 50,920 8.978 8.824 

co | | 
53.147 9.638 9.642 
46.197 | 8.786 | 8.509 

NH, | 

50.862 10.035 9.932 
So, 59.293 9.514 9.528 


(a) Para o CO, e NH, valores de [16]; para o SO, valores de [17]. 
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Oscilador harmônico linear 
Interacção rotação-vibração 
Distorção centrifuga 
Anarmonicidade 


Construções Técnicas ergue no Mondego, esta barragem de abóbodas múltiplas, com BO metros de altura 
e 350 metros de comprimento, que inundará uma área de 2 000 hectares. 


Serão executados 350 000 mº de escavações, utilizados 450 000 mº de betão e 10 000 toneladas de aço. 
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Orientação do projecto de docas secas 


(Estaleiro naval da Setenave, Setúbal) 


RESUMO 


Mencionam-se os actuais condicionamentos relati- 
vos ao projecto de docas secas para construção e para 
reparações navais, com enfase nos imperativos da sua 
integração em estaleiros orientados como verdadeiras 
unidades industriais. Descreve-se sucintamente a orien- 
tação seguida no projecto das docas do novo estaleiro 
português em acabamento no estuário do Sado, Setúbal. 


1 —- Dez ou vinte anos atrás, ainda algumas admi- 
nistrações portuárias consideravam as docas secas como 
«obras públicas», que qualquer um poderia utilizar me- 
diante co pagamento de taxas por uma ocupação avulsa, 
As docas secas cerviam apenas para a conservação & 
reparação dos navios, já que a construção de novos 
cascos ainda se faia em carreiras e não em docas 
atribuídas, exclusivamente ou não, a essa actividade, 
como hojo é prá'ica corrente, pelo menos, para navios 
de grande porte. 

Uma acentuada competividade e uma óntica indus- 
trizl cada vez mais presentes nas actividades da constru- 
ção e da reparação navais, as possibilidades do recurso 
a equipamentos de elevado custo inicial e uma melhor 
compreensão dos benefícios do trabalho organizado con- 
denam hoje tudo O que não seja uma boa integração das 
docas secas em conjuntos sistematizados englobando as 
oficinas e outros meios indispensáveis para as operações 
da um estaleiro completo. As docas deixaram assim de 
ser concebidas com um certo cariz de monumento pe- 
rene, cujos custos e prazos de construção, embora con- 


J. M. Gonçalves Figueira (1) 


SUMMARY 


The main conditions governing shipbuilding and shi- 
prepairing dry docks design is mentioned. Emphasis is 
put on the requirements for integration of these docks in 
shipyard projects conceived as high productive industrial 
units. 

Procedures adopted for the dry docks design at the 
new shipyard at the Setubal estuary (Portugal) are re- 
ported. 


trolados pelos planos e orçamentos de gastos dos di- 
nheiros públicos, não estão sujeitos às duras exigências 
do mercado. 


Os cadernos de encargos a que os promotores da 
construção de novos estaleiros submetem os projectistas 
das docas secas obrigam-nos a considerá-las como um 
equipamento industrial, que como qualquer outro tem de 
ser rentável a tem de se inserir num sistema operacional 
cingindo-se a uma economia de conjunto. Como qual- 
que: outro equipamento, embora habitualmente construi- 
das com materiais de duração praticamente ilimitada, 
o seu custo terá de ser amortizado em curto prazo, pois 
em poucos anos podem tornar-se absoletas por mutações 
técnicas, ou dispensáveis por alterações do mercado 
do transporte maritimo. Os prazos para os estudos e 
construção das docas são muitas vezes drasticamente 
influenciados por condições de oportunidade comercial 
e mutações inesperadas das conjunturas políticas e eco- 
nómicas introduzem factores adicionais de intranquili- 
dade nos que participam na implementação deste tipo 
de obras. (2) 


(1) — Engenheiro civil |. 5. T. consultor na PROFABRIL — Centro de Projectos — Lisboa. 


(2) — No estaleiro da Marqueira, operado pela LISNAVE, por três vezes foi alterada durante a construção a largura das docas 
secas, para atender ao rápido crescimento das dimensões dos navios tanques encomendados a seguir ao encerramento do Canal 
do Suez, Essas alterações foram possíveis porque a concepção estrutural da secção corrente dessas docas — muros laterais em T inver- 
tido independente da laje de fundo —- permitia construir grande parte de um dos muros laterais e grande parte do fundo enquanto se 
aguardava a decisão sobre a posição do segundo muro. 


(3) — Da 1960 até o início da construção do Estaleiro de Setúbal (1972) o número de graneleiros (navios de mais de 10000 tdw 
destinados a graneis lígudios e sólidos) passou de 3042 para 5082, o porte total de 62 X 10º idw para 258 X 10º tdw, o porte médio 
da 20 000 para 45 000 tdw. 
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O crescimento da frota mundial e do porte dos na- 
vios e o abandono da construção em «carreiras», espe- 
cialmenta dos graneleiros (3), tem motivado a partir da 
década de 1960 o estabelecimento de grande número 
da docas de grandes dimensões. Actualmente estão em 
operação, em todo o mundo, cerca de vinte docas para 
navios acima de 300 000 tdw, número que quase dupli- 
cará nos próximos dois anos. 


“4 Sd DA Na ES 4 


...4 as A a A 


A 


a 


túbal), a 45 km a SE de Lisboa destina-se à construção e 
reparação de navios de grande porte, até 700 000 tdw. 
Dispõe de 3 docas secas concebidas nara funcionarem 
como um conjunto coeso, integrado no sistema global 
de produção. Existe uma divisão bem mercada entre 
um sub-sistema de construção de novos navios e im 
sub-sistema de reparação, ocupando, respect vamente, 
a zona W e a zona E do estaleiro. À dispasição adoptada 


+ Fiará dl 
“ 


A gravura mostra a península de Setúbal e os estuários do Tejo e do Sado onde se situam os dois estaleiros portugueses para navios 
da grande porte: LISNAVE, instalado na década de 60 e SETENAVE em acabamento 


A situação que se apresenta vem facultando aos 
projectistas de docas secas um campo de acção vasto 
o pleno da interesse. Mas impõe-lhes atitudes e procedi- 
mentos de profissionalismo, em que engenheiros genera- 
listas e engenheiros especialistas das variadas técnicas 
têm de trabalhar em equipa com os futuros operadores 
do conjunto de instalações em que se integra a obra, 
tirando o melhor partido dos aperfeiçoamentos em cada 
disciplina, mas subordinando a sua participação a uma 
concepção técnica da conjunto e aos imperativos de 
condições de exploração extremamente exigentes quanto 
a produtividade e oportunidade de oferta de serviços. 


2 — O estaleiro que a empresa SETENAVE tem em 
vias de acabamento no estuário do Sado (Porto de Se- 


para o conjunto das docas reflecte o eguilibrio encon- 
trado entre diversos requisitos opostos: especiali:ação 
de funções co integração em cada um dos sub-sistem s, 
procura da máxima utilização de meios comuns, minimi- 
zação de custos da sua própria construção, adaptação e 
um programa de avanço dessa construção favorável ao 
programa de conjunto de todas as obras e equipamentos 
do estaleiro e à lista das primeiras encomendas conse- 
guidas pela empresa promotora, versatilidade de utiliza- 
ção permitindo boas condições de trabalho para várias 
proporções entre as encomendas de trabalho de constru- 
ção de novos navios e de trabalho de reparação. 

Assim, o conjunto das docas situa-se no meio da 
área da produção do estaleiro. Duas delas, em segui- 
mento uma da outra, com o mesmo eixo longitudinal, 
ficam adjacentes à zona das oficinas afectas à constru- 
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1 Doca de construção ZONA DE REPARAÇÃO 
2 Doca de múltipla utilização 10 Caldeiraria ligeira 
3 Doca de reparação 11 Oficinas de mecânica e de tubos 
4 Pontes-cais 12 Escritórios centrais 
5 Cais de descargas do aço 13 Vestiários balneários 
ZONA DE CONSTRUÇÃO 14 Escola 
6 Parque de aço 15 Cantinas 
7 Caldeiraria principal 16 Armazéns gerais 
8 Planos de premontagem 17 Subestação principal 
9 Oficinas de aprestamento e de tubos 18 Central de compressores 


PLANTA GERAL 


O estaleiro ocuna uma área de cerca de 800 000 m? conquistados ao estuário por aterros de areias dragadas. O conjunto das 3 docas 
secas ocupa uma posição central entre a zona W (construção de novos cascos e seu aprestamento) & a zona E (reparação). À zona NE 
à a dos serviços gerais [escritórios centrais, cantina central, serviços médicos, armazém geral, sub-estação prncipal, central de com- 
pressores e central da gases). 
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CORTE TRANSVERSAL DO CONJUNTO DAS DOCAS SECAS 


DC 


O Ra LO 


BUILDING DOCH 


ne 
= 


IR 


ac PR E 


= 
- 
= 
-— 
= 
= 
—— 
mm 


MEM Mir 
rr 


1; Li 


E 


CORTE LONGITUDINAL: Doca de construção e doca de utilização multipl 


cão. À terceira doca, do lado E, confina cem a zona das 
instalações afectas às reparações, 


Das duas docas co-axiais, a de montante estã des- 
tinada a ser a principal zona de montagem final de novos 
cascos: a doca de juzante, além de funcionar como 
eclusa para a comunicação da primeira com o estuário, 
será utilizada quer para montagem, quer para reparação. 
Estas docas constituem uma obra co mcaracterísticas 
muito diferentes do habitual. A estrutura da doca de 
montante é constituída apenas por uma laje encaixada 
entre dois muros, que sobressaem cerca de 6 m do nivel 
terreno circundante. O fundo fica cerca de 2m abaixo 
desse nível, À doca de juzante teria a configuração ha- 
bitual das docas secas se o coroamento dos seus muros 
não tivesse de obedecer ao mesmo nível do coroamento 
da doca adjacente, por forma a possibilitar ao seu en- 
chimento simultâneo e, assim, o seu funcionamento 
como eclusa permitindo transferir para o estuário os 
cascos acabados de construir. O fundo anda pela cota 
— 7 m, cerca de 10 m abaixo do fundo da primeira, per- 
fazendo-se assim uma altura total de 18 m. As dimen- 
sões desta doca, 75 X 450 m, permitem a eclusagem de 
cascos para navios até 700 000 tdw. A doca de mon- 
tante, com a mesma largura, ficou com 420 m de com- 
primento, mas pode ser facilmente ampliada para com- 
primentos praticamente ilimitados. Isto quando interes- 
sar montagem simultânea de navios que requeira tal 
ampliação. 


As duas docas deste grupo são servidas com um 
pórtico para 500 t, cujo vão de 150 m permite varrer o 
terrapleno adjacente do lado W, onde se localizam planos 
de pre-montagem de secções até âquela capacidade de 
elevação. Estes planos recebem por sua vez, transporta- 
das em zorras ou transferidas por guindastes, secções 
do menores dimensões preparadas na caldeiraria prin- 
cipal ou em áreas ao ar livre adjacentes. 


Qualquer das duas docas pode ser compartimentada 
em diversos troços por comportas amovíveis manobra- 
das, em secções, pelo grande pórtico. Uma dessas portas 
à utilizada para fechar o topo posterior da doca de mon- 
tente, quando do alagamento desta doca. A doca de ju- 
zanta dispõe de uma comporta de rebater, utilizada para 
o seu esvaziamento. Quando deva ser cheia a níveis 
superiores aos das águas do estuário, utiliza-se uma 
comporta do tipo das comportas amoviveis referidas, que 
suporta um impulso no sentido de terra para o mar que 
poderá atingir o correspondente a perto de 18 m de água. 


A doca do lado E foi concebida para operações de 
reparação, situando-se na proximidade das respectivas 
oficinas. A sua utilização como doca de construção é 
possivel, mas menos indicada, por obrigar a percursos 
longos na ligação à caldeiraria principal e aos planos 
de pré-montagem a por não estar dotada dos potentes 
meios de movimentação de cargas que guarnecem o 
grupo das duas outras docas. Tem uma configuração tri- 
vial, com o coroamento de nivel com os terrzplenos. 
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A casa das bombas e o centro de comando comuns 
às três docas situam-se junto à margem ,entre as testas 
das docas. À sua capacidade (100 000 m'/h, cerca de 
5000 CV) foi escolhida de acordo com as necessidades 
de rapidez de docagem para serviço de reparação e per- 
mite esvaziar a maior das docas utilizável nesse serviço 
em pouco mais de duas horas. 


A movimentação de cargas nas docas é assegurada 
pelo pórtico de 500t já mencionado e por guindastes 
de lança, parte dos quais circula no terrapleno entre as 
docas e pode servir qualquer delas. Deve sublinhar-se, 
como particularidade da concepção original da doca de 


CORTE DO COROAMENTO DOS MUROS 
DA DOCA DE REPARAÇÃO 


1 ar comprimido 

* água doce 

3 cabos eléctricos 

4 agua, serviço de incêndios 
5 água, serviço geral 
6 água, lastro 
acetilena 

oxigénio 

9 carril de guindaste 
10 cabeço de amarração 
11 quadros eléctricos 


12 caminho da rolamento dos 
carros da alagem de navios 


13 carril dos baileus de pintura 
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construção, a possibilidade de acesso directo ao fundo 
dessa doca de grandes cargas transportadas em veículos 
rodados, 


Além do equipamento de movimentação de cargas, 
quer as operações de construção quer as de reparação, 
exigem um variado conjunto de instalações. A figura em 
que se esquematiza a geometria do coroamento de um 
dos muros das docas e a disposição do equipamento nela 
instalado dá ideia do vulto e da variedade dessas insta- 
lações. Para a movimentação dos navios utiliza-se um 
sistema em que os cabos de alagem puxados por guin- 
chos instalados nos topos das docas são guiados por 
carrinhos que circulam em caminhos de rolamento insta- 
lados ao longo dos muros laterais; sistema que facilita 
o manter o navio bem centrado durante a movimentação 
(o que é muito vantajoso, em especial, quando a folga 
aos muros laterais é apenas de alguns decimetros). Às 
docas são servidas por redes de energia preparadas para 
as necessidades das operações fabris, da manobra do 
equipamento específico das docas e do abastecimento 
dos navios em reparação (estes exigindo correntes com 
diversas características, consoante as normas adoptadas 
pelos seus armadores). A água é fornecida por várias 
redes conforme o seu destino, destacando-se, como valor 
do caudal disponivel, a rede de água de lastro, dimen- 
sionada para o enchimento rápido dos navios tanques em 
ensaio. As operações fabris, além do sistema eléctrico, 
são apoiadas pelas redes de ar comprimido e de gases. 
As docas são ainda dotadas de caminhos de rolamen- 
to para um baileu destinado a trabalhos ligeiros no casco 
dos navios (pinturas, etc.) sem necessidade de escadas 
ou andaimes que embaraçassem a circulação no fundo 
das docas. 


A elevada parcela com que a estrutura das docas 
concorre para o custo total de um estaleiro de constru- 
ção/reparação navais, justifica grande apuro na investi- 
gação das condições naturais do local e na escolha dos 
tipos de estrutura e dos métodos do seu dimensiona- 
mentc. Do vulto destas obras dão ideia os 180 000 m' 
de betão, 14 000 t de varão de aço, 150 000 m” de su- 
perfícies de moldes, e as 6500 estacas (12m X 7) 52 em, 
moldadas no terreno) aplicados no grupo das três docas 
em apreço, 


As docas estão engastadas num maciço de areias 
limpas, muito permeáveis (coef. de perm. 3 X 10(2m/s) 
constituído por uma camada natural aluvial, de espes- 
sura praticamente ilimitada pelo que respeita às condi- 
cões impostas às obras, completada por um aterro de 
cerca de 56 m executado com areias dragadas nessa 
mesma camada. As características desse maciço foram 
cuidadosamente investigadas: sondagens com amostra- 
gem e SPTs, análises granulométricas, sondagens com 
penetrómetros estáticos e dinâmicos, medições directas 
da permeabilidade e ensaios à escala natural das estacas 
de ancoragem. Durante a obra os aterros foram controla- 
dos, quer quanto às compacidades obtidas, quer quanto 
à eventual inclusão de lentilhas lodosas, 


A alta permeabilidade do terreno constitui uma 
forte determinante na escolha dos tipos de estrutura, 
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MANOBRA DE LANÇAMENTO DE UM NOVO CASCO 


i-— na doca de construção está um casco concluído a uma secção de rá em construção; a doca-eclusa está livre; ambas as docas 


vazias. 


ii — as docas estão alagadas até ao nivel máximo; o novo casco está sendo transferido para a doca-eclusa; aproveita-se o alagamento 


para mudar a posição da secção de rê. 


rt — a água é descida ate deixar seca a doca de construção, acertando com o nivel do estuário na doca-eclusa. 


iv -— aberta a porta, o casco à rebocado para fora, deixando a doca-eclusa livre para qualquer utilização. 


quo suportam pressões hidrostáticas praticamente cor- 
respondentes ao nivel da maré no estuário, e dos proce- 
dimentos de construção. Para as duas docas mais pro- 
fundas, a construção ao abrigo dum importante sistema 
de rebaixamento do nivel freático impunha-se como 
solução não contestável e a ancoragem dos fundos mos- 
trou-se francamente mais vantajosa do que o contraba- 
lançar as pressões hidrostáticas exclusivamente com 
o peso das estruturas. Com a solução adoptada para a 
doca de construção, fundada a uma cota onde as impul- 
sões hidráulicas são insignificantes, foram possíveis uma 
importante economia no custo da estrutura e a execução 
em prazo muito curto. Vantagens estas que, juntas às 
que o esquema oferecia à operação do estaleiro, mencio- 
nadas atrás, foram confrontadas favoravelmente com 
gastos de instalação e inconvenientes de exploração de 
sinal contrário, nomeadamente, acréscimo de custo da 
doca-eclusa (muros laterais sobre-elevados, porta suple- 
menta suportando impulsos de dentro para fora). Assim, 
a doca de construção pode ser iniciada na fase inicial do 
rebaixamento do nivel freático da área ocupada pelo 
conjunto das três docas (operação que ocupou 6 dos 24 
meses que levaram a concluir as estruturas do conjunto); 
pôdo ser adoptada uma fundação directa sobre os ater- 
ros compactados, executada com um rebaixamento das 
dguas praticamente equivalente apenas à amplitude da 
maré, solução que além da poupança do custo das esta- 


cas de ancoragem, a que se recorreu para as restantes 
docas, permitiu uma redução suplementar de prazo. 


O apuro do reconhecimento geotecnico consentiu 
dimensionar as estruturas com muita confiança no valor 
das solicitações do meio encaixante assim como permitiu 
dimensionar com justeza as obras provisórias (enseca- 
deira a sistema do rebaixamento de águas). Escolhidas 
as estacas de betão moldadas no terreno como meio de 
ancoragem contra as sub-pressões criou-se um problema 
de dimensionamento assás complexo e trabalhoso na me- 
dida em que os requisitos de minimização de custos im- 
punham um dimensionamento à justa. À situação batia 
nos seguintes pontos: estacas especialmente escolhidas 
com vista a funcionarem como tirantes verticais, mas 
também solicitadas (alternadamente) a compressões sig- 
nificativas, grupos de estacas inclinadas para ancoragem 
dos esforços devidos à impulsão (de dentro para fora) 
na grande porta da doca-eclusa; grande diversidade de 
combinações de solicitações, nomeadamente no caso 
de ensaios de enchimento dos navios tanques, A utiliza- 
ção do computador (IBM 1130, com traçador de grá- 
ficos) deu uma saída satisfatória para as dificuldades, 
permitindo abordar os cálculos da laja da fundo como 
lajo apoiada simultaneamente em dois meios elásticos 
com módulos de elasticidade diferentes (terreno e sis- 
tema de estacas), considerar as variações de secção da 
lajo o a sua integração (para a secção ccrrente) numa 
peça em U englobando os dois muros laterais. 
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O empreendimento do estaleiro de Setúbal foi lan- 
cado em princípios de 1971, realizando-se então os 
estudos técnicos preliminares e a escolha do terreno, 
numa área do estuário do Sado a conquistar por aterro. 
Às dragagens e aterros ocuparam o ano de 1972 e du- 
rante este ano desenvolveu-se a maior parte do trabalho 
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de projecto. A construção das docas iniciou-se com o ano 
da 1973. Em princípios de 1974 já a doca de constru- 
ção permitia, embora não terminada, o início da monta- 
gem do primeiro casco. Em Janeiro de 1975 foram dadas 
por terminadas as duas docas mais baixas, iniciando-se 
a remoção da respectiva ensecadeira. 


SECÇÃO TRANSVERSAL DA ESTRUTURA DA DOCA DE CONSTRUÇÃO 


Fundada directamente sobre o terreno arenoso, a secção transversal compreenda duas peças laterais, separadas pela viga de quilha 
central. Cada peça lateral é moldada com 30 m de comprimento. Nas restantes duas docas a secção transversal corrente é um U mo- 
nolítico, fundado no conjunto terreno + estacas, moldado em troços também de 30 m 


TÉCNICA 430 


87 


Feiras, Conferências, Exposições e Congressos 


OBJECTIVOS E PROGRAMA 


DIDACTA — European Educational Materials Fair 


Ill Conference ON ELECTRODE PHENOMENA 
IN GAS DISCHARGES 


V Simpósio lberoamericano de Catálise 


me =— ES mto dE TESE ea o add 


International Trade Fair 


First International Timplate 
Conference 


| 


Lisboa 


O o 


LOCAL, DATA, ORGANIZAÇÕES 
E INFORMAÇÕES 


Basilca 
23 a 27 de Março - 76 


Lisboa 
Maio - 76 


E E = - E — 


25 a 30 de Julho - 76 
Fundação Gulbenkian 


Elmia, Jônkóping, Sweden 
27 de Setembro a 1 de Outubro - 76 


Londres 
5 a 8 de Outubro - 76 
In Research Institute 


E E E — 4-0 == TS = — e 


Madrid 
11, 12 e 13 de Novembro - 76 


é FUNDAÇÕES ESPECIAIS 
. ESTUDOS GEOTÉCNICOS 


: CAPTAÇÕES DE ÁGUA 


Witt R.DO ARSENAL, 146 -2º 
tino Telf. 320208 -364010-360437 


Cio LISBOA -2 


FABRICANTE DE ELÉCTRODOS E ARAME | 
PARA SOLDADURA AUTOMÁTICA 
EQUIPAMENTO PARA SOLDADURA  ELEC- 
TRICA E OXI-ACETILÊENICA 
ACESSÓRIOS 
ad 

MATERIAL AFROVADO POR 

AMERICAN BUREAU OF SHIPPING 
BUREAU VERITAS 
DET NORSKE VERITAS 
GERMANISCHER LLOYD 
INSTITUTO DE SOLDADURA 

LLOYD'S REGISTER OF SHIPPING 


* 


SEDE, FÁBRICA E LABORATÓRIOS 
Rua Vice-Almirante Azevedo Coutinho, 4 
VENDA NOVA — AMADORA 


Telefs. 9776321 .971733- 976471 - 970613 | 


DELEGAÇÃO NO NORTE 
Rua Delfim Ferreira, 525 
PORTO 
Telef. 621 24 


poTalia bel 


2 * VOLUME 


Pelo Engenheiro 


À. 0. AERCI 


Preço de cada volume 


180$00 


Desconto: 10º/, aos assinantes 


Pedidos à Técnica 


TECNICA X 


convectores 
eléctricos 


EFAGEL 


calor de prestígio 


dia lie 


TÉCNICA IX 


Resumo dos Artigos publicados na «Técnica» n.º 430 


Ano L —- Novembro 1975 


C. D. U.621-135:532.5 
ANTONIO A. 5. DE BARROS 


ANALISE TEÓRICA DAS CARACTERISTICAS MECANICAS 
DE UM ESCOAMENTO RADIAL NUMA TURBOMAQUINA 
CENTRIFUGA 


Técnica N.º 430 — L— 11-1975, p. 47 a 66 


Neste trabalho são apresentadas, de uma forma condensada, 
as deduções das equações gerais das velocidades e acele- 
rações a energia num escoamento sem perdas no interior 
da roda de uma turbomáquina radial, É, seguidamente, 
feita a aplicação das equações estabelecidas a vários es- 
coamentos: escoamento a energia constante e escoamentos 
com aumento de energia (velocidade relativa constante, 
inclinação radial constante e ná plana) e, finalmente, escoa- 
mentos com entrada e saida impostas: pá circular e variação 
linear da inclinação tangencial. 


C. D. U. 536:533.27.5 


J. &. MARTINHO SIMÕES 


a 
C. A. NIETO DE CASTRO 


TERMODINAMICA MOLECULAR DO GAS PERFEITO 
Técnica N.º 430 — L— 111975, p. 67 a 80 


Apresenta-se uma revisão crítica da determinação teórica 
de propriedades termodinâmicas de gases perfeitos. Com 
baso na Termodinâmica Estatística e em dados fornecidos 
pela analise espectroscópica indicam-se métodos de cálculo 
da várias contribuições para os valores das funções termo- 
dinâmicas. Finalmente, aplica-se o estudo às moléculas de 
CO. NH, e SO, e comparam-se resultados calculados de 
S e Cp com valores experimentais. 
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Cc. D. U. 627.31/32 
J. M. GONÇALVES FIGUEIRA 


URIENTAÇÃO DO PROJECTO DE DOCAS SECAS 
[ESTALEIRO NAVAL DA SETENAVE, SETUBAL) 


Tócnica N.º 430 — L-— 114975, p. 81 a 88 


Mencionam-se os actuais condicionamentos relativos ao 
projecto de docas secas para construção e para reparações 
navais, com enfase nos imperativos da sua integração em 
estaleiros orientados como verdadeiras unidades industriais. 
Descreve-se sucintamente a orientação seguida no projecto 
das docas do novo estaleiro português em acabamento no 
estuário do Sado, Setúbal. 
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Synopsis of articles published in «Técnica» n.º 430 


U. D. C. 627.31/32 
J. M. GONÇALVES FIGUEIRA 


REQUIREMENTS FOR DRY DOCKS DESIGN 
(SETENAVE SHIPYARD, SETUBAL) 


Tecnica N.º 4390 — L— 1141-1975, p. B1 a 88 


The main conditions governing shipbuilding and shiprepai- 
ring dry docks design is mentioned. Emphasis is put on 
tha requirements for integration of these docks in shipyard 
projects conceived as high productiva industrial units. 
Procedures adopted for rhe dry docks design at the new 
shipyard at the Setubal estuary (Portugal) are reported. 


L - November 1975 


U. D. C. 621135:532.5 
ANTONIO A. 5. DE BARROS 


THEORETICAL ANALYSIS OF THE MECHANICAL FEATURES 
OF THE FLOW IN A CENTRIFUGAL TURBOMACHINE 


Técnica N.º 430 — L=— 11-1975, p. 47 a 66 


Tris study establish the velouty and aceleration equations 
for the flow in a centrifugal turbomachine, wich are applied 
to some types of flow: without and with energy increase. 


U. D. C. 536:533.27.5 
J. à. MARTINHO SIMÕES 
e 
C. A. NIETO DE CASTRO 


MOLECULAR THERMODYNAMICS 
OF THE PERFECT GAS 


Técnica N.º 430 — L— 11-1975, p. 67 a 80 


à critical review about the theoretical determination of 
theromdynamic properties of perfect gases is made. Based 
in Statistical thermodynamics and in data derived from 
spectroscopic analysis, methods of calculus of several con- 
tributions for the values of the thermodynamic functions 
are indicated. Finally, the study is applied to the calcula- 
tion of moar entropy and heat capacity of simpla po- 
Iyatomic molecules, viz CO, NH, and S0,. À comparison 
with experimental values is made. 
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